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RESUME 


Bien que les véhicules électriques soient de plus en plus présents sur les routes et que la 
majorité des manufacturiers automobiles préparent leur stratégie d’électrification, un élé- 
ment en particulier représente toujours un défi de taille. Malgré les avancées des dernières 
années sur la technologie des accumulateurs au lithium-ion, produire un bloc-batterie 
contenant assez d'énergie et de puissance à un coût raisonnable pour rivaliser avec les 
moteurs à combustion n’est pas une tâche simple. De plus, avec le nombre de modèles, de 
format et de chimie qui ne cesse d'augmenter, évaluer toutes les possibilités devient une 
táche titanesque sans l'aide d'un ordinateur. Traditionnellement, l'essai-erreur était trés 
populaire pour le design de bloc-batterie, mais gráce aux puissances de calcul disponible 
aujourd'hui, l'optimisation devient une solution de plus en plus populaire pour améliorer 
chaque aspect de celui-ci. On remarque que la premiere phase du design, potentiellement 
la plus importante, oü le choix des cellules et de configurations séries-paralléles est fait 
est souvent négligée. En automatisant le calcul des solutions faisables et l'optimisation 
pour obtenir le maximum d'information concernant le design, beaucoup de temps peut 
étre sauvé avant de passé aux analyses d'ingénierie approfondie. Ce projet de recherche 
propose une nouvelle méthodologie de dimensionnement de bloc-batterie voulant répondre 
à ce probléme. 


La méthodologie présentée dans ce projet de recherche est divisée en 2 parties. Tout 
d'abord, un Constraint Satisfaction Problem (CSP) est utilisé en combinaison avec la re- 
présentation N,/N, pour produire l'ensemble des solutions faisables selon les contraintes 
énoncées par le concepteur. Cette premiére étape sert principalement à éliminer tous les 
designs qui ne respectent pas les contraintes du véhicule. Par la suite, le résultat de ce 
CSP agit comme intrant pour un Constrained Optimization Problem (COP) qui ajoute 
une couche d'optimisation pour extraire le maximum d'informations sur les meilleurs de- 
signs potentiels. Pour obtenir les solutions non dominées, le probléme multi-objectif est 
transformé en mono-objectif par des fonctions scalaires. En combinant plusieurs sommes 
pondérées évaluées avec des poids différents ainsi que l'évaluation de métriques L4 et Loo, 
cela permet d'obtenir un front de Pareto plus diversifié. De plus, les métriques sont calcu- 
lées à partir d'un point de référence qui permet de refléter les préférences du concepteur 
dans l'algorithme d'optimisation. Pour faciliter le travail du concepteur, la méthodolo- 
gie est implémentée à l'intérieur d'une interface graphique programmée avec MATLAB. 
Celle-ci permet de contróler directement le processus de sélection des contraintes et l'op- 
timisation en plus de donner accés à sept modeles de cellules différentes. Pour démontrer 
l'efficacité de la méthodologie, celle-ci est appliquée sur un cas réel de dimensionnement 
pour une moto électrique de course. Chaque étape y est décrite pour mieux comprendre le 
processus d'obtention de l'ensemble des solutions faisables et d'optimisation. Les résultats 
démontrent que la sélection d'une configuration séries-paralléles n'est pas triviale. Dans 
le front de Pareto obtenu, la métrique L4 a été choisi pour sa nature compensatoire, pour 
un design équilibré de bloc-batterie pour une moto électrique. 


Mots-clés : Multi-objectif, blocs-batterie, véhicule électrique, optimisation, CSP, COP 


ABSTRACT 


Battery pack sizing based on series/parallel configurations as a constrained 
optimization problem 


While electric vehicles are more and more present on the roads and the majority of automo- 
tive manufacturers are preparing their electrification strategies, one element in particular 
still represents a major challenge. Despite major advancements in the past years in the 
lithium-ion battery technology, producing a battery pack containing enough energy and 
power with a reasonable cost to rivalize with internal combustion engines is not a triv- 
ial endeavour. Moreover, with the ever increasing number of cell models, format and 
chemistry, evaluating every possible combination is becoming increasingly difficult to do 
without automation. Traditionally, trial and error was the method of choice for the design 
of battery packs but with the help of the increasing computing power available today, 
optimization is becoming a very popular method to better every aspect of the design. Its 
been noted that the in the first phase of the design, potentially the most important, the 
choice of cell and series-parallel configuration is often neglected. By automating this part 
of the process and by giving as much information as possible to the designer, a lot of 
time can be saved before going into in-depth engineering analyses. This research project 
proposes a new methodology for the sizing of battery packs. 


The methodology presented in this research project is divided into two parts. First, a 
Constraint Satisfaction Problem (CSP) is used in combination with the N,/N, represen- 
tation to produce the feasible set of solutions that meets the selected constraints. This 
first step is used to eliminate all the designs that do not meet the vehicles requirements. 
Thereafter, the results of the CSP is used as inputs for a Constrained Optimization Prob- 
lem (COP) that adds an optimization layer to extract as much information on the best 
available designs. To obtain the non-dominated solutions, the multi-objective problem is 
scalarized to a mono-objective problem. By combining many weighted sums evaluated 
with different sets of weights and by also evaluating the Lı and La metrics, this allows 
to obtain a diversified and broad Pareto front. Moreover, the metrics are calculated in 
respect to a reference point defined with engineering values to reflect the preferences of 
the designer in the optimization. To ease the work of the designer, the methodology is 
implemented within a Graphical User Interface programmed with MATLAB. This allows 
direct control over the selection of constraints and the optimization process while allow- 
ing the choice of seven different cells. To demonstrate to efficiency of the methodology, 
it is applied to a real case of battery pack sizing for an electric motorcycle. Every step 
is detailed to understand how the algorithm works, to understand how the feasible set 
solutions are obtained and how they are optimized. The results shows that the selecton 
of the series-parallel configurations is not a trivial choice. The choice was made based on 
the Lı metric for its compensatory nature to obtain a balanced design. 


Keywords: Multi-objective, Battery pack, Electric vehicles, optimization, CSP, COP 
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CHAPITRE 1 
INTRODUCTION 


1.1 Mise en contexte et problématique 


Depuis déjà quelques années, le marché des véhicules électriques prend de plus en plus 
d’ampleur. Bien que ces véhicules existent depuis plus d’un siècle, l’évolution de la techno- 
logie des batteries au cours des dernières années a permis d’augmenter drastiquement leurs 
performances. Le regain d'intérêt envers cette technologie n'est pas un hasard. Les chan- 
gements climatiques et la déplétion des ressources de pétrole forcent le passage vers des 
énergies renouvelables. L'électrification des transports se présente comme étant la solution 


la plus prometteuse pour la diminution des gaz à effet de serre (GES). 


Cependant, une électrification compléte des transports ainsi qu'une transition vers une 
production d'énergie intégralement renouvelable est encore loin d'étre réalisée. Cette tran- 
sition est limitée par un facteur commun aux deux secteurs, soit la batterie. Malgré les 
avancées des dernières années, cette technologie n'est pas capable de rivaliser avec son 
compétiteur, le pétrole. Plusieurs laboratoires tentent de perfectionner la chimie des cel- 
lules pour en améliorer la puissance, l'énergie et la sécurité. D'autres se concentrent plus à 
l'intégration de ces cellules dans des batteries pour la création de systémes d'accumulation 
d'énergie (SAE). Il s'agit d'un combat sur plusieurs fronts pour amener cette technologie 


à un niveau pouvant dépasser les énergies fossiles. 


La conception d'un SAE est un processus d'ingénierie complexe, onéreux et qui requiert 
énormément de temps. En alliant génie mécanique, électrique, chimique et informatique, 
il s'agit d'un effort multidisciplinaire pour pouvoir tenir compte de tous les phénoménes 
concurrents se produisant dans la batterie. De plus, une multitude de paramétres doivent 
faire partie de l'étude de faisabilité pour obtenir un produit répondant aux besoins du 
client. Le chemin pour arriver à une batterie optimale est donc trés difficile et comporte 
plusieurs défis. Pour ce genre de probléme, les ingénieurs ont de plus en plus recours à 
l'ingénierie assistée par ordinateur (CAE). L'accessibilité de la puissance de calcul permet 
aujourd'hui d'effectuer des simulations et des opérations trés complexes beaucoup plus 
facilement qu'il y a quelques années. Cette technologie offre un niveau de fidélité et de 
précision incomparable aux méthodes traditionnelles. C'est pourquoi un outil informa- 


tique d'aide à la conception semble étre une contribution pour résoudre le probléme. En 
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réduisant les délais de livraison et en facilitant le processus, beaucoup d’argent pourra être 
épargné et de nouveaux concepts pourront être itérés plus rapidement. De plus, en utilisant 
des techniques d'optimisation, il sera plus simple pour l'ingénieur de trouver le compromis 
optimal pour les besoins de la situation. La validation reste une partie du processus qui 
est difficile et complexe. Mais en l'intégrant à un outil informatique, le concepteur pourra 
directement comparer et valider ses concepts pour s'assurer d'avoir le meilleur des com- 
promis. Donc non seulement ceux-ci se feront plus rapidement, mais l'ingénieur sera en 


mesure de livrer un meilleur produit. 


1.2 Description du projet de recherche 


Aujourd'hui avec la démocratisation des véhicules électriques, les manufacturiers se lancent 
dans la course pour obtenir le plus de parts du marché possible. Il en découle une vaste 
sélection de cellules lithium-ion de différente chimie, de différent format, performance, prix 
et disponibilité. Ceci entraine l'impossibilité pour un concepteur de tenir compte de tout ce 
qui se retrouve sur le marché sans l'aide d'un logiciel automatisant une partie du processus. 
Le méme probléme est rencontré lorsqu'il doit évaluer ensuite les requis pour le véhicule 
et tenir compte des différents compromis possibles selon les options disponibles. Que ce 
soit d'évaluer les méthodes de refroidissement, la structure du bloc-batterie, l'arrangement 
des cellules, la täche est maintenant trop complexe pour étre effectué efficacement avec 
les méthodologies traditionnelles. De plus, un élément souvent oublié est la maniére dont 
les différentes configurations du bloc-batterie sont évaluées de façon qualitative. Pourquoi 
devrait-on prendre un arrangement de cellule plutót qu'un autre ? Comment celui-ci peut-il 


étre optimisé ? 
On peut donc en déduire 2 questions de recherche : 


1. Est-il possible d'élaborer une méthodologie pour améliorer, faciliter et accélérer le 


processus d'ingénierie du dimensionnement du bloc-batterie ? 


2. Comment peut-on quantifier la qualité des résultats obtenus avec la nouvelle métho- 


dologie en vue d'une comparaison avec les approches traditionnelles ? 


De ces questions découlent les objectifs et contributions originales de ce mémoire défini 
comme suit. Tout d'abord, le développement d'un outil de dimensionnement permettant 
de comparer de multiples cellules entre elles et en considérant le Phase Change Composite 
(PCC) comme matériel de refroidissement. Il automatise une partie du processus, à l'aide 
d'une interface graphique, l'évaluation des contraintes choisies et l'évaluation des objectifs 


du concepteur pour obtenir un front de Pareto selon ses préférences. Ce qui laisse plus 
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de temps et de liberté au concepteur pour itérer son concept et étudier les différentes 
possibilités. Par la suite, la qualité des solutions est évaluée à l’aide d’une méthode d’op- 
timisation pour donner au concepteur une vue d’ensemble sur les compromis possibles en 


plus d’évaluer lesquels seront susceptibles d’être plus performants selon les requis. 


1.3 Plan du mémoire 


Le mémoire est divisé en 5 chapitres et une annexe de la manière suivante. Le chapitre 1 
présente une introduction et une mise en contexte par rapport au bloc-batterie lithium-ion 
pour véhicule électrique, les défis actuels de l’industrie automobile ainsi que les méthodes 
de dimensionnement de ceux-ci. Également, les objectifs de ce projet de recherche sont 
énoncés et les contributions originales sont mises en valeur. Le chapitre 2 passe en revue 
l'état de l'art sur les blocs-batterie lithium-ion, les différentes méthodologies de dimension- 
nement ainsi que l’optimisation multi-objectif. Le chapitre 3 présente une méthode ma- 
thématique utilisant la théorie des ensembles et le problème de satisfaction de contraintes 
pour calculer les configurations possibles d’un bloc-batterie selon les requis clients. En 
utilisant la représentation série-parallèle, permettant de voir l’ensemble des solutions fai- 
sables pour la conception d’un bloc-batterie. Le chapitre 4 construit sur le méthodologie 
présenté au chapitre 3 en apportant un cadre d’optimisation et le développement d’un 
logiciel personnalisé pour le dimensionnement de bloc-batterie. En étendant le problème 
de satisfaction de contraintes vers un problème d'optimisation contraint, l'ensemble des 
solutions faisables est optimisé selon les préférences du concepteur. Ceci s'effectue en com- 
binant des sommes pondérées de fonctions scalaires avec le calcul des métriques L4 et Læ 
(distance Manhattan et distance Chebyshev respectivemment). Ces distances sont évaluées 
par rapport à un point de référence calculé à partir de valeur d'ingénierie permettant de 
refléter les préférences du concepteur. Le chapitre 5 discute des principales conclusions 
du mémoire et propose des pistes intéressantes pour des travaux futurs. Finalement, l'an- 
nexe présente la description pour l'ajout d'une cellule à l'interface graphique ainsi que la 


maniére de sauvegarder et de charger une configuration. 


CHAPITRE 1. INTRODUCTION 


CHAPITRE 2 


Etat de l'art 


2.1 Systeme d’accumulation d'énergie 


On définit un système d’accumulation d’energie (SAE) comme étant le processus de conver- 
sion de l'énergie électrique provenant d'un réseau électrique dans une forme ou elle peut 
être stockée pour être reconvertie lorsque nécessaire [7]. Ces systèmes possèdent plusieurs 
applications tant dans les appareils portables, dans le transport ainsi que pour les res- 
sources d'énergies renouvelables. Ces dernières sont en forte croissance et les SAE se dé- 
marquent comme solution idéale pouvant stocker l'énergie étant donnée la nature impré- 
visible de production d'électricité provenant de celles-ci. C'est pourquoi les applications 
stationnaires de bloc-batterie seront un marché très intéressant d’ici quelques années. On 
peut classifier les SAE dans quatre catégories, soit électriques, mécaniques, chimiques ou 
thermiques [8]. Pour les applications véhiculaires, deux catégories se démarquent, soit les 
systèmes d’accumulation d’énergie électrique comme les super-condensateurs et ceux de 
type électrochimique plus communément appelés batteries. Cette recherche se concentre 


principalement sur les accumulateurs électrochimiques de type lithium-ion. 


Les accumulateurs électrochimiques ont pour fonction de transformer l'énergie produite 
par des réactions électrochimiques en énergie électrique. Lors d’une décharge, des réactions 
se produisent aux deux électrodes, soit l’anode et la cathode, ce qui produit un flux 
d'électrons vers le circuit externe et décharge le générateur. Ces réactions peuvent être 
réversibles ou non. Lorsque le processus n’est pas réversible, on parle de pile, comme 
les piles alcalines. Elles ne peuvent donc pas être rechargées. Lorsque le processus est 
réversible, on parle d’accumulateurs, comme les cellules lithium-ion. En appliquant une 
source d'énergie extérieure, on peut recharger la batterie [7]. Depuis déjà plusieurs années, 
les accumulateurs se présentent comme étant le SAE de prédilection pour l'électrification 


des transports. 


2.1.1 Batterie lithium-ion 


Un bloc-batterie est l'assemblage d'un ensemble d'accumulateurs connectés entre eux pour 
former une source de tension et de capacité qui répond aux besoins d'une application. Pour 


un véhicule électrique, la batterie est généralement divisée en trois niveaux, soit la cellule, 
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le module et la batterie. La cellule représente la plus petite unité de la batterie, il s’agit 
d’un accumulateur électrochimique possédant la plus basse tension selon sa chimie [9]. Les 
cellules sont ensuite interconnectées en série et en parallèle pour produire la configuration 
répondant aux besoins de l’application. Chaque module possède son propre système de 
gestion. Les contraintes matérielles vont engendrer une limitation du nombre de cellules 
en séries. Les modules sont ensuite interconnectés en série et en parallèle pour former 
la batterie. Généralement, la batterie sera protégée par une structure de plastique ou de 
métal comprenant un système de gestion de la batterie complète ainsi qu’un système de 


refroidissement [10]. 


La modélisation d’accumulateurs au lithium, telle que la cellule, permet d’étendre le mo- 
dèle vers une batterie virtuelle [11]. Pour que ce modèle soit fidèle à la réalité et pour être 
en mesure de l’étendre vers un bloc-batterie virtuel, plusieurs facteurs sont à considérer. 
La prochaine section examine les composants qui forment la batterie et leurs impacts sur 


la conception et le dimensionnement. 
Format des cellules lithium-ion 


Il existe plusieurs formats de cellules pour répondre aux besoins des différents marchés. 
Trois formats se distinguent de ce lot : les cellules en pochettes, les cellules prismatiques et 
les cellules cylindriques. Les cellules en pochettes sont fabriquées en empilant des plaques 
de matériaux actifs à l’intérieur d’un emballage non rigide et flexible. Celles-ci requièrent 
une attention particulière lors de manutention et ont besoin d’une enveloppe mécanique 
rigide supplémentaire. Quant aux cellules prismatiques et cylindriques, elles sont consti- 
tuées d’un boîtier robuste qui protège les matériaux actifs et qui empêche l’expansion de 


la cellule tout au long de sa durée de vie [10]. 


La cellule étant l’unité de base d’une batterie, elle dictera la direction de la conception 
de la batterie. Il est primordial de bien étudier les besoins pour choisir le format adéquat. 
Les présents travaux permettront une évaluation de l’impact des différents formats sur la 


conception de la batterie. 
Configuration des cellules 


Le format de la cellule à un impact sur la configuration de celles-ci à l’intérieur de la 
batterie. Les cellules en pochettes et prismatiques sont généralement placées de façon 
linéaire alors que les cellules de type cylindriques peuvent aussi être disposées de façon 
quinconce comme illustré à la figure 2.1. La topologie choisie ainsi que l'arrangement série 


et parallèle aura un impact non seulement sur la méthode de refroidissement, mais aussi 
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Figure 2.1 Arrangement quinconce 


sur le vieillissement de la batterie [12]. La difference de capacité des cellules introduite 
lors de la fabrication réduit la capacité totale de la batterie. Pour réduire la dispersion 
statistique de cette valeur et pour maximiser la capacité de la batterie, une configuration de 
connexion en série de modules parallèle est à privilégier. De plus, la façon dont les cellules 
sont placées dans la batterie, selon leur différence de caractéristiques, possède aussi une 
influence sur la capacité de la batterie [13]. L’arrangement et la topologie des cellules sont 
des points cruciaux de la conception d’une batterie. Pour optimiser la configuration séries- 
parallele, il est primordial de bien considérer ce facteur étant donné qu'il peut affecter les 


performances, la gestion thermique, le vieillissement et la capacité de la batterie. 
Systéme de gestion de la batterie 


Le systéme de gestion de la batterie est un composant essentiel pour le bon fonctionnement 
de la batterie ainsi que pour la sécurité. Il protége les cellules de surtension, de sous- 
tension, des demandes de courant trop élevées lors de recharge ou décharge ainsi que des 
fluctuations de températures trop grandes. Son autre fonction principale est le balancement 
des cellules. Deux techniques sont présentées dans la littérature soit le balancement passif 
(ou dissipatif) et l'actif (ou non dissipatif). Lors d'un balancement passif, un regroupement 
de résistances est utilisé pour retirer l'énergie des cellules possédant un plus haut état de 
charge pour égaliser les tensions avec les cellules ayant le bas état de charge. Il ne s'agit 
pas d'une méthode trés efficace énergétiquement et elle génére de la chaleur. Pour le 
balancement actif, l'énergie est transférée des cellules ayant un état de charge élevé vers 


les cellules ayant un état de charge bas, à l'aide de condensateur ou d’inductance [10]. 
Méthodes de refroidissement 


Pour obtenir des performances optimales lors de l'utilisation d'une batterie dans un véhi- 


cule, un certain intervalle de température doit étre respecté. Il est essentiel de conserver 
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les cellules lithium-ion autour sous les 60 °C pour éviter les dommages, améliorer la durée 
de vie de la batterie ainsi que conserver ses performances. Pour respecter ces contraintes, 
la température des cellules devrait être maintenue autour de 25 °C [14]. Les différentes 
méthodes de refroidissement ont toutes leurs forces, mais aussi leurs faiblesses. Il est donc 
primordial d'avoir la technique de refroidissement adéquate pour évacuer la chaleur gé- 
nérée lors de l’utilisation du véhicule, mais aussi conserver la température uniforme. Le 
choix de la méthode influence la conception de la batterie et est important à considérer 


pour obtenir une solution optimale au problème. 


L’approche la plus simple est le refroidissement à l’air. En forçant un écoulement d’air dans 
la batterie, on peut limiter les hausses de température. Par contre, un usage très agressif 
de la batterie ou dans des conditions d’opérations du système à température élevée, ce 
genre de refroidissement peut causer de larges zones non uniformes de température qui 
pourraient engendrer des problèmes. La méthode de refroidissement par liquide est plus 
efficace, mais complexe, coûteuse, et possède des risques de fuites. Une approche par 
matériaux à changement de phase peut permettre de limiter la température maximale 
atteinte par la batterie, mais peut ajouter beaucoup de masse et volume [2]. Le choix 
de la méthode de refroidissement a un impact considérable sur l’aspect thermique de la 
batterie, mais aussi sur le volume, la masse et le coût. Il faut donc en tenir compte lors de 


la conception pour choisir la méthode optimale à la situation. 
Dimensionnement d'une batterie 


L'intégration d'une batterie dans un véhicule électrique représente un défi considérable. 
Il faut tout d'abord s'assurer de bien définir les limites d'utilisation, les conditions envi- 
ronnementales dans lesquelles la batterie devra opérer ainsi que les exigences techniques. 


Pour ce faire, il faut bien étudier les différentes caractéristiques principales d'un SAE soit : 


— ]a tension — le volume 

— l'énergie — le poids 

— ]a puissance — le coût 

— la temperature — la durée de vie 


Certaines de ces caractéristiques agiront comme filtres, puisqu’elles sont fixes selon les 
besoins du client comme la tension par exemple. Ceci permet de réduire dès le départ le 
nombre de solutions. D’autres seront beaucoup plus complexes à déterminer. Elles pour- 
raient même comporter plusieurs solutions possibles. L’optimisation de ces caractéristiques 


peut être très difficile étant donné que plusieurs des objectifs sont contradictoires. Un bon 
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exemple, si l’on optimise la dur&e de vie, la conception sera probablement beaucoup plus 
coûteux et si l’on optimise le coût il aura une durée de vie réduite. Il n’est pas trivial de 
trouver l’arrangement optimal d’une batterie, on recherche plutôt le compromis optimal. 
Bien que l’expérience d’un ingénieur peut s’avérer très utile dans ce genre de situation, 
l'utilisation d'un processus d'optimisation numérique multi-objectif devient indispensable 
lorsque l’on veut explorer tout l’espace des solutions. La majorité des méthodologies de 
conception de bloc-batterie utilisent des méthodes d'optimisation. Cependant, il n'existe 
pas une méthode structurée infaillible à toutes les situations, mais plutôt une multitude 
de méthodes pour répondre à la majorité des problèmes. La prochaine section décrit l’état 
de l’art en ce qui concerne ces méthodes et surligne les opportunités d'améliorations de 


chacune d’entre elles. 


2.1.2 Méthodes de dimensionnement de batterie lithium-ion 


Une fois les besoins analysés et le cahier des charges établi, il faut choisir une méthodologie 
appropriée pour la réalisation du projet et obtenir un produit final répondant aux besoins 
énoncés. Étant donnée la non-trivialité de la conception de bloc-batterie, cette étape est 
cruciale. Traditionnellement, l’approche était de modéliser par conception assisté par or- 
dinateur le bloc-batterie et à l’aide d’essais et erreurs on pouvait trouver un arrangement 
satisfaisant les requis. Cependant, si celui-ci ne répond pas aux requis suite aux essais, le 
processus devenait très lourd sur le plan du temps [4]. C’est pourquoi plusieurs méthodes 
ont été proposées pour assister, automatiser ou améliorer le processus de conception des 
blocs-batterie [1, 2, 3, 4, 5, 15, 16, 17, 18]. Ces méthodes peuvent être regroupées en quatre 


grandes catégories : 

1. Modélisation mathématique 

2. L’optimisation 

З. Représentation N,;/N, 

4. Logiciels de conception 
Cette sous-section explore chacune des grandes catégories pour mieux comprendre les 
forces mais aussi les possibilités d'améliorations. 


Méthodologies traditionnelles basées sur la modélisation 


La trés grande majorité des méthodologies de conception de bloc-batterie lithium-ion vont 
utiliser un modéle mathématique de batterie. Les non-linéarités présentes augmentent 
grandement la complexité, presque jusqu'à l'impossibilité de traduire ces comportements 


en équations mathématiques. C'est à ce moment que les modéles mathématiques de bat- 
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terie entrent en jeux. Une méthode de conception qui minimise le coût total et tente de 
maximiser la durée de vie a été proposée par [1]. Celui-ci utilise un modèle électrochimique 
nommé Enchanced Single Particle Model (ESPM). Le modèle est utilisé pour effectuer des 
simulations de différentes batteries à différentes températures. La figure 2.2 résume le 


processus de décision pour la conception du bloc-batterie. 


Active 
temperature 
control available? 


Yes 


Determine the expected battery 

operation temperature and use it 
with figure 3(c) to select the 

optimal pack size for maximum 


Maximum life [ec ume RS Minimum cost 


cost ($/day) or 
maximum life 
(vears) 


life and minimum cost 


Battery pack temperature should 
be 25° Celsius. Use it with figure 
3(c) to select the optimal pack 
size for maximum life 


E cost of battery pack using figure 3(c)($) 
Tes life of battery pack using figure 3(c) (days) Er 


Cr is cost ($/day) to maintain the battery temperature. Get 
the Т. for overall minimum cost. Use Т with figure 3(c) 


to select the optimal pack size for minimum cost 


Figure 2.2 Processus de conception pour la sélection de la batterie optimale 
tiré de [1] 


L'utilisation d'un modéle électrochimique seul ne permet pas de tenir en compte tous 
les aspects d'une batterie. De plus, il est souvent trop précis et requiert beaucoup de 
travail pour bien le paramétrer. Bien que cette méthode fournisse une bonne idée de ce 
que la batterie pourrait avoir comme volume et comportement thermique, bien d'autres 
paramétres restent à étre étudiés pour obtenir un concept optimal. Il est cependant un 


des rares à considérer le vieillissement dans sa méthodologie. 


Il est aussi possible de modéliser une cellule mathématiquement à l'aide de modéle équi- 
valent électrique, tel que [16]. Le modéle empirique permet de simplifier la modélisation 
et de réduire l'effort requis pour la paramétrisation. Cette méthodologie se concentre à 
trouver des lois sans-dimensions de mise à l'échelle pour les caractéristiques d'une batterie. 


L'utilisation d'une base de données contenant les différents paramétres du modéle dans 
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différentes situations obtenus expérimentalement simplifie la paramétrisation du modele et 
permet une modélisation plus fiable. La méthode est divisée en deux parties, tout d’abord 
il faut trouver la meilleure cellule et ensuite effectuer le dimensionnement de la batterie. 
La deuxième étant basée sur la première. En fournissant des spécifications à l’algorithme, 
celui-ci tentera de trouver une cellule qui permettra d’obtenir les cibles. Il s’agit d’une 
méthode simple, mais efficace. Cependant, celle-ci amène une surestimation du nombre de 
cellules et donc de la masse. Bien que l’on puisse simuler le comportement de la tension, 
cette méthode ne prend pas en compte l’aspect thermique et le vieillissement de la batterie. 


L’aspect économique n’est pas considéré dans cette méthode. 


Les méthodologies basées sur des modèles sont très fréquemment utilisées. Elles permettent 
d'obtenir généralement de bonnes estimations, mais ne sont pas complètes. La lourdeur du 
calcul imposée par la simulation limite la précision des modèles. En effet, un modèle que 
l’on qualifierait de parfait serait beaucoup trop coûteux à simuler sur les ressources infor- 
matiques d’aujourd’hui. L’exploration complète de l’espace de solutions serait beaucoup 
trop longue. C’est pour cela que la majorité d’entre eux vont se concentrer sur un aspect 
soit électrochimique ou électrique pour faire un modèle. Il s’agit souvent de la méthode 
de base puisqu'elle est commune à la majorité des autres méthodologies. Elles vont prin- 
cipalement ajouter des outils pour améliorer la méthodologie et ainsi obtenir de meilleurs 


résultats. Par exemple, avec des techniques d'optimisation. 
Dimensionnement basé sur la configuration séries-parallèle 


De nouvelles manières de représenter l’arrangement d’un bloc-batterie de façon graphique 
sont de plus en plus populaires. Par exemple, la représentation série-parallèle, ou №, /№,, 
qui affiche les résultats en utilisant le nombre de cellules en série et en parallele [19, 20, 21]. 
Cette représentation peut étre trés utile pour bien illustrer de maniére concise les diffé- 
rentes possibilités pour visualiser les différents concepts répondant aux requis. En utilisant 
cette représentation, on peut afficher par exemple l'ensemble des solutions faisables selon 
les requis d'un véhicule |20]. De plus, cette méthode s'adapte bien à l'intégration d'un 
algorithme d'optimisation pour passer d'un ensemble de solutions faisables à un ensemble 


de solutions non dominées pouvant être examinées selon les compromis à faire. 
Méthodologies utilisant l'optimisation 


Les ressources informatiques d'aujourd'hui permettent d'effectuer des calculs trés com- 
plexes. Pour cette raison, les ingénieurs se tournent de plus en plus vers l'optimisation 
numérique. Étant donné les multitudes de solutions possibles, l'objectif est de trouver le 


meilleur compromis. C'est ce que plusieurs proposent d'utiliser pour optimiser la concep- 
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tion des blocs-batterie lithium-ion. La méthode proposée par [2] utilise l'algorithme Multi- 
objective Particule Swarm Optimization (MOPSO) pour optimiser la conception du sys- 
téme de refroidissement en concordance avec la réponse thermique de la batterie. On peut 
voir le processus à la figure 2.3. Le processus est semblable à ce que l'on retrouve à la 
section précédente, mais apporte des améliorations. Cette méthode permet d'explorer plus 
de solutions en obtenant le front de Pareto. De cette fagon, l'ingénieur peut choisir plus 
facilement le meilleur compromis. Cette méthodologie se concentre sur la modélisation 
thermique de la batterie, donc certains aspects sont laissés de cóté. De plus, l'utilisation 
du logiciel COMSOL Multiphysics(R) pour la co-simulation augmente les coüts d'utilisation 


et la complexité. 


LiBESS Optimal Design 


Framework 
-— n 
Step 1: Cell 
Cell Selection & =» Eee 
N Pack Sizing 


Application 
Requirements 


d Step 2: 
Technical Constraints Arrangement Pattern & 


М В SR, 
Cooling System Solution J Models | 
& 


Electrical Use Pattern V 
“a ^ Optimization 
Step 3: Tools 
Environmental Arrangement & Cooling System 
Conditions 


x Optimization Variables P 


Step 4: 
Long Term Evaluation & 
V Actual Cost Value Estimation J 


Figure 2.3 Processus de conception de bloc-batterie optimisé tiré de [2] 


On peut retrouver aussi des méthodes utilisant des méthodes d'optimisation hybrides qui 
combine optimisation locale et globale pour profiter des avantages des deux approches. 
C'est le cas de la méthode proposée par [3] qui propose une méthode pour minimiser la 
masse, le volume et le coüt d'un bloc-batterie avec une méthode hybride. Une premiere 
optimisation globale est effectuée pour trouver les variables entiéres optimales. Celles-ci 
sont ensuite utilisées comme point de départ pour l'optimisation locale qui a pour but 
de localiser précisément l'optimum global. On peut voir le processus à la figure 2.4. On 
remarque que cette méthode permet de prendre en compte une multitude de facteurs 
durant l'optimisation pour réduire la masse, le volume et le coüt en plus de pouvoir 
effectuer une simulation de la batterie gräce aux modeles. Il s'agit d'un outil qui permet 
d'obtenir une trés bonne idée des différents concepts dans les phases initiales de conception. 


Par contre, il s'agit d'une méthode trés complexe et difficile à reproduire. En plus, l'article 
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stipule qu’il faut une semaine complete de calcul pour arriver à l’optimum. Le vieillissement 


n’est pas inclus dans cette méthode. 
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Figure 2.4 Processus d'optimisation hybride de batterie lithium-ion tiré de [3] 


Les méthodologies utilisant l’optimisation apportent une couche de plus que celles utilisant 
seulement un modèle. On peut arriver plus rapidement à un concept optimal, mais ce 
sont des processus beaucoup plus complexes et peu conviviaux. De plus, ces méthodes 
ne permettent pas beaucoup de flexibilité à l’ingénieur. П peut difficilement faire des 
changements de variables et comparer les concepts. Les études de sensibilités sont donc 
très longues à exécuter avec ces processus. Avec ces problèmes en tête, certains ont décidé 


d'aller vers l'implémentation d'outils informatique pour se simplifier la vie. 
Méthodologies basées sur des outils informatiques 


La virtualisation des prototypes est un phénoméne trés fréquent aujourd'hui en ingénie- 
rie. On peut travailler sur un concept sans méme avoir besoin d'un prototype matériel. 
On épargne ainsi temps et argent. Le concept de bloc-batterie n'y échappe pas. Il existe 
plusieurs logiciels d'ingénierie assistée par ordinateur qui permettent de calculer le com- 
portement électrique et électrochimique des batteries. Cependant, les coüts d'installations 
sont trés élevés et l'utilisation de ces logiciels requiert du personnel formé [4]. Ces logiciels 


sont souvent trés spécialisés et ne traitent qu'une partie du processus de conception. C'est 
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pourquoi plusieurs chercheurs tentent d’offrir une solution qui épaulera l’ingénieur tout au 


long du processus grâce à des outils informatiques. 


Certains proposent d'utiliser des outils et des méthodes de prototypage virtuel pour évaluer 
tôt dans le processus le comportement thermique de la batterie [4]. L'objectif est de réduire 
le temps de développement et les coûts. La plateforme conçue est illustrée à la figure 2.5. 
À l’aide d'une interface graphique, l'utilisateur peut configurer le type de test à effectuer 
sur la batterie. Il peut choisir le type et la chimie de la cellule. Il peut aussi configurer la 
manière dont les cellules sont placées dans la batterie et choisir le type de refroidissement 


utiliser. 


/ BatteryPack \ l Electric layout 
>! 
H " Specification ) K input | Cell datasheet 


Testing : i| Cell Electrical 
Configuration Г тезї 


HD Tools : 


m ip simplified battery 
Pack: CFD р model and 
Thermal Analysis : cooling analysis 


Figure 2.5 Plateforme de conception de bloc-batterie lithium-ion tiré de [4] 


Bien que l’interface graphique permette de faciliter grandement le processus de conception, 
celui-ci n’est pas complet. Le coüt et le vieillissement ne sont pas pris en compte lors 
des analyses, elles sont principalement concentrées sur le comportement thermique. De 
plus, l’outil n’est pas complètement indépendant, il a besoin d’un logiciel commercial 
d’analyse de la dynamique des fluides. L’architecture de cette plateforme vient compliquer 
le processus de conception, possède un coût de démarrage élevé et requiert un employé 


formé à opérer le logiciel. 


Une automatisation du processus de conception électrique et mécanique a été proposée 
par [5]. À l'aide d'une interface graphique, voir la figure 2.6, il permet au concepteur de 
générer un prototype mécanique virtuel qui rencontre des spécifications comme la puis- 
sance, l'énergie ainsi que la facilité d'assemblage. Par contre, plusieurs simplifications sont 
effectuées et la batterie n’est pas étudiée en profondeur. On obtient principalement une 
bonne idée de la configuration mécanique ainsi que sa facilité d'assemblage. Le comporte- 
ment thermique et électrique durant son utilisation et lors de recharge et décharge n'est 


pas considéré. 
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Figure 2.6 Exemple d’interface graphique d’outil de dimensionnement tir& de 


[5] 


L’intégration des méthodes de conception soit par modéles ou optimisation dans des outils 
informatiques est un pas dans la bonne direction. La simplification du processus d’ingénie- 
rie est nécessaire pour réduire les coüts et le temps de développement. Cependant, aucune 
méthode jusqu’à maintenant ne permet de soutenir l’ingénieur tout au long du processus. 
La combinaison de l’optimisation, de la modélisation et du génie logiciel semble être l’ave- 
nue la plus prometteuse. En prenant les forces de chacune des méthodes actuelles et en 
tentant de corriger les faiblesses, l’objectif est de développer une meilleure méthodologie. 
Pour cela, il faut explorer chacun des trois aspects en profondeur et trouver la meilleure 


combinaison. 
Logiciels commerciaux de dimensionnement 


Avec les grandes compagnies automobiles qui développent maintenant des stratégies d’élec- 
trification ou de production de masse de véhicules électriques, de plus en plus de compa- 
gnies offrent des solutions logiciels pour le dimensionnement complet de bloc-batterie. Que 
ce soit pour la conception mécanique, l’optimisation ou la simulation ces logiciels sont de 
plus en plus complets, mais aussi complexes. On cite principalement les compagnies comme 
ANSYS, Gamma Technologies et COMSOL Multiphysics. Cette solution comme méthode 
de dimensionnement est cependant trés coüteuse avec des licences annuelles onéreuses ainsi 


que le requis d’avoir du personnel spécialement formé pour l’utiliser. 
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2.2 Programmation par contraintes 


Lorsque l’on veut utiliser un logiciel pour résoudre un problème de conception, il faut 
pouvoir traduire le problème en quelque chose de compréhensible par l’ordinateur. Dans 
le cas de la conception d’un bloc-batterie, le défi est principalement de comparer toutes 
les possibilités disponibles comme les cellules, l'arrangement et les différentes méthodes de 
refroidissement. Par contre, effectuer cette comparaison à la main est pratiquement impos- 
sible, dû au nombre en croissance constante de nouveaux modèles de cellules, mais aussi 
des combinaisons quasi infinies possibles dans l’arrangement de celles-ci. Une méthode 
particulièrement efficace pour résoudre des problèmes combinatoires complexes est la pro- 
grammation par contraintes. Cette méthode consiste à utiliser les liens logiques entre les 
variables, leurs domaines et les contraintes les reliant pour définir un ensemble de solutions. 


Deux méthodes en particulier sont présentées dans cette section. 


2.2.1 Problème de satisfaction de contraintes 


L'utilisation du probléme de satisfaction de contraintes (CSP) pour des applications en 
ingénierie regagne en popularité dans les dernières années [22, 23]. Apparue vers la fin 
des années 1980, cette méthode de résolutions de problèmes combinatoires est issue de la 
programmation logique et de l'intelligence artificielle [24]. Dans un CSP, un ensemble de 
variables possède un domaine fini et est relié entre elles par un ensemble de contraintes. 
Plusieurs types de CSP existent, soit le MAX-CSP, où on tente de satisfaire un nombre 
maximal de contraintes dures en transgressant le moins de contraintes molles possible [25]. 
Une autre définition serait le problème de satisfaction de contraintes pesées (WCSP) où 
l'on tente de minimiser le poids total des contraintes non satisfaites [26]. Le but principal 
est donc d'effectuer une recherche qui à la fin ne donnera que les solutions possibles. Pour 
qu'une solution soit valide, elle doit étre compléte et constante : elle doit donc respecter 


le domaine fini de la variable ainsi que toutes les contraintes associées respectivement. 


2.2.2 Probléme d'optimisation contraint 


Bien que d'obtenir un ensemble de solutions faisables soit trés intéressant, l'ajout d'une 
couche d'optimisisation s'avére étre important pour améliorer les concepts en permettant 
de mieux comprendre les compromis pour respecter les requis. De par sa similarité de 
structure, le CSP peut étre étendue à un probléme d'optimisation contraint (COP) [16]. Au 
lieu de considérer seulement une seule fonction objectif comme dans un CSP, un ensemble 
de plusieurs fonctions objectifs est utilisés pour la résolution du COP. Cela permet de 


prendre le résultat initial de solutions faisables pour ensuite trouver les meilleures solutions 
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disponibles. Bien qu'il soit possible de faire des analyses d'optimisation exhaustive avec 
des logiciels trés puissants, il est parfois mieux d'obtenir une premiére bonne solution 
plus rapidement que l'optimum, mais qui est non-dominée [27]. De plus, il est difficile de 
bien comprendre les compromis à faire et leurs impacts dés le départ. En favorisant les 
itérations gráce à une optimisation plus rapide, le concepteur peut tester différentes options 
rapidement et diriger la recherche vers un résultat plus souhaitable selon le cahier des 
charges. Cette méthode peut étre utilisée avec n'importe quel algorithme d'optimisation. 
Le choix de celui-ci devient donc une étape importante pour une résolution efficace qui 


donne un front de Pareto avec des solutions non dominées diversifiées. 


2.3 Optimisation multi-objectif 


Il est indéniable que les systémes d'ingénierie sont de plus en plus complexes et donc de 
plus en plus difficiles à concevoir et produire. De plus, avec l'utilisation des logiciels de 
conception assistée par ordinateur, on numérise fréquemment les produits à développer. 
Par la suite, on veut optimiser la solution pour délivrer le meilleur produit possible. C'est 


pourquoi les ingénieurs se tournent vers le domaine de l'optimisation. 


L'optimisation est un processus par lequel on veut répondre à la question «s'agit-il de 
la meilleure option possible ?», contrairement à une approche traditionnelle en ingénierie 
oü la question est plutót «est-ce suffisant ?», pour un probléme donné. Il est important 
de bien définir le terme «meilleur», car il est défini par des critéres établis à priori selon 
le contexte du probléme et son application. Le choix de ces critéres est primordial pour 
obtenir la solution optimale. On doit tout d'abord comprendre quels sont les paramétres 
les plus intéressants à faire varier. Il s'agit des variables du probléme. Ensuite, il faut 
déterminer quelles sont les limites dans lesquelles on peut les faire varier. On définit l'espace 
de recherche pertinent. Pour concevoir les fonctions objectifs, il faut se demander quels 
sont les objectifs et comment on peut les exprimer sous forme mathématique. Finalement, 
il faut décider de la méthode d'optimisation à utiliser. L'avantage principal de ce processus 
est de pouvoir générer rapidement une multitude de concepts différents et de les analyser 


par rapport à nos critéres définis a priori. 


De plus en plus, on doit considérer plusieurs objectifs contradictoires simultanément. De 
surcroít, avec les nouvelles avancées en intelligence artificielle ainsi que la montée crois- 
sante des puissances de calcul, une multitude de nouvelles possibilités sont disponibles en 
optimisation multi-objectif. Le couplage de ces méthodes avec des modéles virtuels est de 


plus en plus fréquent avec une adoption croissante des méthodes d'ingénierie prédictive. 
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Encore une fois, la connaissance de son problème ainsi qu'une formulation claire et précise 
est primordiale pour une optimisation efficace. Diverses méthodes sont proposées dans la 
littérature et tous ont une particularité ou un type de problème particulier auquel elles 
peuvent répondre. Il n’existe pas de méthode parfaite pour toutes les situations, mais bien 
une multitude de méthodes pour répondre à tous les types de problèmes. Étant donnée la 


nature de la recherche, seules les méthodes d'optimisation numérique sont considérées. 


2.3.1 Méthodes locales 


Les méthodes basées sur le calcul local du gradient se divisent en deux parties, soit les 
méthodes itératives et les méthodes directes. Ces méthodes ont pour but de converger 
vers le minimum local le plus rapidement possible. Une méthode itérative populaire est 
l'algorithme de descente raide (steepest descent) qui minimise une fonction réelle et diffé- 
rentiable. On procède par améliorations successives en calculant la dérivée de la fonction 
objectif pour se déplacer vers le point stationnaire de façon proportionnelle à la pente. 
L'itération est représentée de la façon suivante où 7,4, est le prochain point à calculer, ть 


le point actuel, o le pas de l’itération courante et па la dérivée de la fonction objectif : 


df 
Teri = Ir — Ор 2.1 
k+1 k KT (2.1) 
Cette méthode n’est pas utilisable dans un contexte où le problème n’est pas différentiable. 
Pour l'utiliser dans le contexte de cette recherche, on doit absolument le coupler à une 


autre méthode. Il en va de méme pour les autres méthodes dites continues et itératives. 


La deuxiéme catégorie de méthodes d'optimisation locale est celle des méthodes directes. 
La méthode simplex qui a été développé par [28] permet l'optimisation d'une fonction 
objectif méme si elle contient des discontinuités ou qu'elles ne sont pas différentiables. On 
cherche un point stationnaire voisin directement à partir de points évalués de la fonction 
objectif à des endroits stratégiques. Si la fonction objectif contient n variables, on évalue 
n+1 points équidistants pour former un tétraédre à n dimensions. Une fois l'évaluation des 
points effectués, le point le plus éloigné par rapport au point stationnaire est réfléchi pour 
construire le prochain tétraédre. Cette méthode est aussi connue comme la programmation 


linéaire. 
2.3.2 Méthodes globales 


Les méthodes globales (ou aléatoires) se divisent elles aussi en deux catégories soit les 


méthodes de marche aléatoire et les méthodes de recherche stochastiques. Ce sont des 
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méthodes plus flexibles que les méthodes locales. Elles permettent d'éviter de rester pris au 
piége dans un minimum ou un maximum local. La marche aléatoire est un processus simple, 
mais peu efficace. On calcule la fonction objectif pour un grand nombre de points choisi de 
facon aléatoire dans l'espace de la fonction et on retient la meilleure solution. Les méthodes 
dites stochastiques utilisent un processus aléatoire pour générer des perturbations pour 
permettre d'obtenir plusieurs points de maximum ou minimum locaux. On utilise ces 
méthodes lorsque la fonction objectif est non convexe, non différentiable et non continue. 
Ce sont des méthodes simples à mettre en oeuvre, par contre elles ont un coüt de calcul 
trés élevé. Il faut donc les ressources informatiques en conséquence puisque pour chaque 


solution, une simulation et une évaluation complete de la fonction de coüt sont nécessaires. 
Algorithme génétique 


L'une des méthodes globales populaires aujourd'hui est l'algorithme génétique. Elle se base 
la théorie de l'évolution des espéces et la théorie de la génétique pour parcourir l'espace 
objectif [29, 30, 31]. Les individus, soit les solutions possibles à un probléme posé repré- 
sentent la population. Celle-ci évolue de fagon itérative par les générations jusqu'à ce que le 
critére d'arrét soit atteint. Lors de chaque génération, plusieurs opérateurs sont appliqués 
aux individus pour engendrer la nouvelle population. Tout d'abord, un croisement d'indi- 
vidus sera effectué pour engendrer de nouvelles solutions, soit les enfants. Un opérateur 
de mutation peut étre appliqué pour introduire un bruit et sortir d'un potentiel minimum 
local. Son róle est de modifier un individu de fagon aléatoire pour en produire un nouveau. 
Ensuite, une évaluation des enfants est faite selon une fonction de performance pour déci- 
der quelles solutions seront conservées. On termine avec le remplacement des solutions les 
plus faibles par de meilleures solutions. Un des algorithmes génétiques les plus utilisés est 
celui proposé par |29], nommé NSGA-II. Sa stratégie élitiste, un triage de solutions non 
dominées rapide et une façon efficace de gérer les contraintes en ont fait un algorithme 
de choix. Bien que cette méthode soit efficace, lorsque le nombre de variables est élevé et 
que les simulations pour chaque individu sont longues, le temps de calcul explose. N'ayant 
pas de mémoire, il est possible qu'un individu soit calculé plus d'une fois. Par contre, sa 


simplicité est un atout lors de l'implémentation. 
Recuit simulé 


La méthode du recuit simulé développé par [32] se base sur le phénomène physique portant 
le méme nom. Principe tiré de la métallurgie, il s'agit de la tendance des atomes à s'or- 
ganiser naturellement dans une configuration qui minimise l'énergie lorsque le matériau 


fondu est refroidi lentement. On associe la fonction objectif à l'énergie du systéme, donc 


20 CHAPITRE 2. ETAT DE L’ART 


une solution proche du minimum global aura une faible énergie. Cet algorithme se base sur 
deux résultats de physique statistique. Premièrement, lorsqu'un système atteint l'équilibre 
thermodynamique à une température donnée T, la probabilité de celui-ci de posséder une 
énergie donnée Е, est proportionnelle au facteur de Boltzmann où kg désigne la constante 


de Boltzmann : 


d Biker (2.2) 


Deuxièmement, cet algorithme utilise le critère de Metropolis [33] qui effectue une trans- 
formation élémentaire sur le système et si cette transformation permet de diminuer la 
fonction objectif, alors on accepte la transformation. Dans le cas où celle-ci ferait aug- 
menter la fonction objectif, elle serait rejetée. On applique ce critère de façon itérative 
pour obtenir une séquence de configurations qui représente une chaîne de Markov. Ce 
formalisme permet de dire qu'avec une chaîne de longueur infinie, le système atteint un 
équilibre thermodynamique à la température considérée. Une fois l’équilibre atteint, on 
abaisse quelque peu la température et on refait la chaine de Markov pour ce palier. On 
répète le processus jusqu'à ce que le système devienne immobile, autrement dit, qu'il ny 


a plus de diminution de l’énergie. 
Recherche par tabous 


La recherche par tabous proposée par [34] est un outil trés puissant pour l’optimisation 
de problème combinatoire. Cette méthode utilise une opération appelée mouvement pour 
définir le voisinage d’une solution donnée. On peut le voir comme étant une méthode 
itérative qui, en faisant des mouvements répétitifs d’une solution à un autre, peut explorer 
l'espace solution. À la manière de la descente du gradient. On atteint l'objectif en utilisant 
une liste de grandeur finie de mouvements interdits, appelés tabous, dérivée de l'historique 


récent de la recherche. 


Les deux composants principaux de la recherche par tabous sont la liste de restriction 
et le critére d'aspiration de la solution associé à ces restrictions. La liste de restrictions 
tabou représente la mémoire adaptative puisque quelques attributs sont fixes tant qu'ils 
sont dans la liste. Les listes enregistrent les mouvements dans l'ordre dans lequel ils sont 
faits. Si un nouvel attribut entre dans la liste, le plus vieux est retiré. Le choix de la 
longueur de la liste est critique à l'algorithme. Le critére d'aspiration peut passer outre 
les restrictions. Si un mouvement donné est interdit, mais que le critére d'aspiration est 


satisfait, alors ce mouvement peut étre réactivé. Le critére le plus utilisé est celui qui 
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enlève un tabou lorsque le mouvement essayé donne une meilleure solution que la meilleure 
solution obtenue jusqu'à maintenant. Cette méthode pourrait être très intéressante dans 
le cas de cette recherche. La mémoire permet d'améliorer l'efficacité, mais étant donnée la 
nature du probléme, le concept de voisinage se préte trés bien au jeu. En effet, les concepts 
étant trés similaires, une exploration de ce type pourrait s'avérer plus judicieuse qu'une 


exploration complétement aléatoire comme l'algorithme génétique. 
Essaim de particules 


La méthode d'optimisation par essaim de particules multi-objectif a été proposée par 
[35] comme étant une extension de la méthode proposée par [36] qui est inspirée de la 
chorégraphie d'une volée d'oiseaux. Cette méthode utilise une population dans laquelle le 
comportement de chaque individu est affecté par le meilleur individu local, dans un certain 
rayon, ou bien le meilleur individu global. Cette performance est évaluée de facon similaire 
à un algorithme évolutionnaire et on peut faire la correspondance entre l'ajustement des 
individus à l'opérateur de croisement. La différence majeure est que cet algorithme utilise 
le potentiel de vol à travers l'espace de solution pour en accélérer la convergence. De plus, 
cette méthode garde en mémoire l'expérience de vol de chacune des particules pour pouvoir 
obtenir les solutions non dominées optimales. L'attrait principal de cette méthode est sa 


rapidité de convergence. 


2.3.3 Méthodes hybrides 


Les méthodes locales et globales possédent leurs forces et leurs faiblesses. Certains cher- 
cheurs tentent de mitiger ceux-ci en combinant les deux pour améliorer la convergence, 
mais aussi pour éliminer certains problémes inhérents ont ces méthodes [3]. On retrouve 
de plus en plus de couplage de méthode d'optimisation pour tenter de mitiger les faiblesses 
individuelles |37]. Le couplage reste une nouvelle facon de procéder qui n'est pas encore 
trés utilisée due à sa complexité et le temps de calcul ajouté par rapport à l'utilisation 


d'une seule des deux méthodes. 


2.3.4 Méthodes bi-niveaux 


Certains problèmes d’optimisation possèdent des comportements intrinsèques difficiles à 
écrire en équation. Par exemple, l’optimisation du concept d’un bloc-batterie tout en 
maximisant sa durée de vie. Intégrer le vieillissement dans l’équation n’est pas trivial. 
Une solution intéressante à ce problème est l’optimisation bi-niveaux [38]. Cette méthode 


consiste en deux problèmes d’optimisation emboités l’un dans l’autre. On obtient donc une 
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fonction objectif dans la boucle supérieure et une autre fonction objectif dans la boucle 
inférieure. Ce genre d’approche est de plus en plus utilisée pour régler des problèmes comme 
la tarification de l'électricité [39] et méme le concept de bloc-batterie [6]. La méthode 


proposée par |6] est présentée à la figure 2.7. 


Application Scenario 


ESS and Vehicle Driving cycle 
Parameters profile 


Cost Function Constraints 
Minimization of installation + ESS and Vehicle Models, 
running economical costs. Vehicle power demanded. 


Minimum Batteries + 

Maximum SCs Solution 
Find minimum» j^ with 

Ns = ng using DP 


Battery Only Solution 
Find minimumn ;^with 
пус = 0 using DP 


Outer Loop 


+ 


Define admissible search space 
Battery and SC combinations, between Battery Only and 
Minumim Batteries + Maximum SCs scenarios 


No. of Cells Combinations 


For each combination of 
batteries and SCs: 
Apply Dynamic Programming 
Use battery’s static and dynamic, linear and non-linear 
models and SC’s linear static model. 
From the computed 
solution obtain: 


Optimal Solution 
Battery and SC combination with best economical benefit 


Inner Loop 


d4- — — — — 4-—-—- — — — — — — — — + 


Energy Management 


(minimum installation + running costs) 


Figure 2.7 Algorithme d’optimisation bi niveaux tiré de [6] 
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2.3.5 Approches scalaires 


Bien qu'il est très intéressant de tenir compte de tous les objectifs en même temps, cela est 
souvent un défi de taille. C’est pourquoi, dans certains cas, il est préférable de transformer 
un problème d'optimisation multi-objectif en un probléme mono-objectif à l’aide d'une 
méthode scalaire. Le front de Pareto peut alors être calculé par des méthodes comme la 
somme pondérée ou bien la minimisation d’une distance à un point de référence choisit 


par le concepteur [40]. 
Somme pondérée 


Il existe une multitude d’approches scalaires et la méthode la plus souvent utilisée est 
certainement la somme pondérée. Principalement pour son intuitivité et sa facilité d’im- 
plémentation [41, 42]. Le principe de la somme pondérée est d'appliquer des poids différents 
à chacune des fonctions objectifs et ensuite faire la somme de ceux-ci pour obtenir une 
valeur scalaire. Ceci permet d'évaluer chacune des solutions non dominées sur la même 
base. Pour conserver la condition d'optimalité de Pareto, la somme totale des poids doit 
toujours étre égal à 1. De plus, pour obtenir le pointage agrégé de toutes les solutions, les 
sommes pondérées seront minimisées, c'est-à-dire que si une fonction doit être maximisée, 
le coefficient sera multiplié par -1 |43]. La définition de la somme pondérée est donnée 
à l'équation (2.3). Il existe aussi plusieurs maniéres de définir comment les poids seront 
attribués [44], mais lorsque l'on veut étudier un front de Pareto plus large, l'utilisation de 
valeurs étendues est préférable. Un probléme souvent rencontré avec ce genre d'approche 
est la différence d'unités ou d'ordre de grandeur entre les différentes variables. En utili- 
sation une méthode de normalisation, on peut s'assurer d'obtenir une comparaison plus 
égale [45]. 


Epp Sou Uo: > фу =] (2.3) 


Distances L, et Ly 


Une autre manière d’utiliser une approche scalaire à un problème d’optimisation mulit- 
objectifs est d’utiliser le concept de distance. En définissant un point de référence que le 
concepteur aimerait atteindre en utilisant des valeurs d’ingénierie, il est possible d’inté- 
grer à l'algorithme d’optimisation les préférences de celui-ci. Ce qui permet d'obtenir des 
solutions non dominées plus prés de ce dont le concepteur à besoin. Les deux métriques 


les plus présentent dans la littérature sont la L4 (distance Manhattan) et la Læ (distance 
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Chebyshev). On utilise la distance Manhattan principalement pour compenser les diffe- 
rences absolues entre chacune des fonctions objectives par rapport au point de reference 
pour obtenir une solution plus balancée [46]. La distance Chebyshev est quant à elle non 
compensatoire puisqu'elle ne considère que la plus grande difference entre les objectifs et 
le point de référence. On considére la plus grande déviation, donc le pire scénario |47]. 
En combinant ces deux métriques, il est possible d'ajouter des solutions non dominées 
différentes au front de Pareto et ainsi obtenir plus d'information quant au probléme à 


résoudre. 


Il existe pratiquement autant de méthodes d'optimisation de problémes à optimiser. Choi- 
sir la plus appropriée est donc une opération complexe. Les méthodes globales sont pro- 
metteuses, puisque l'on veut principalement dégrossir le probléme et orienter l'ingénieur 
dans la bonne direction. Cependant, ces méthodes possédent généralement de longs temps 
de calcul. C'est pourquoi une approche scalaire semble étre la meilleure option. De par la 
simplicité de l'implémentation, leur flexibilité dans le nombre d'objectifs possible ainsi que 
leur temps rapide d'exécution, cette méthode semble étre la plus appropriée pour l'appli- 
cation de ce projet de recherche. De plus, il sera plus facile d'incorporer les préférences 
du concepteur à l'aide de cette méthode. Une fois que le concepteur a optimisé son pro- 
bléme et sélectionné un concept, il doit s'assurer que celui-ci effectuera le travail avant de 
construire un prototype matériel. La prochaine section fait l'état de l'art des méthodes de 


simulation et de validation de la conception de bloc-batterie lithium-ion. 


2.4 Simulation et validation du concept 


Une étape importante dans le processus de conception de bloc-batterie est bien sur la 
validation. Il faut pouvoir démontrer que le concept retenu rempli les requit ou bien étre 
capable de comparer plusieurs concepts pour en ressortir le meilleur. L'idéal est de conce- 
voir une batterie virtuelle pour permettre la simulation du cycle qui lui sera appliqué dans 


la réalité pour obtenir un portrait réaliste de son comportement potentiel. 


2.4.1 Modeéles de batterie lithium-ion 


Le domaine de la modélisation de batterie est vaste. Le modèle à utiliser dépend du niveau 
de détails nécessaire énoncé dans la formulation du probléme. Un modele peut étre de trés 
haut niveau pour simplement valider la configuration mécanique de la batterie ou bien trés 
précise pour explorer les interactions électrochimiques lors d'une recharge. Le résultat final 


dépendra toujours de la fiabilité et de la qualité du modéle. Plus un modele sera précis, 
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plus il sera complexe et demandera plus de ressources. Il faut donc trouver l'équilibre selon 
le probléme à résoudre. Il existe plusieurs types de modèles de batterie dans la littérature 


qui permettent d'étudier l'information à différents niveaux. 


Les modèles stochastiques [48] et analytiques [49, 50] bien que parmi les plus simples, ne 
permettent pas de modéliser les non-linéarités et de tenir compte de l'aspect thermique 
et du vieillissement. Il existe des modéles modélisant un aspect particulier de la batterie 
comme les modèles thermiques [51], les équivalents électriques [52] ou bien électrochimiques 
[53]. Bien qu'ils soient trés efficaces dans leur domaine respectif, ils sont souvent trop 
complexes pour étre intégré dans une méthode de dimensionnement ou ne sont pas assez 
complets pour tenir compte de tous les aspects [54]. Les cellules lithium-ion regroupent 
plusieurs phénoménes physiques concurrents, soit électriques, thermiques et chimiques. Il 
est donc impossible de modéliser une batterie à partir d'un seul type de modele, il faut 
utiliser une combinaison de ceux-ci. On remarque que les entrées de certains modéles 
sont les sorties de d'autres. Il est donc possible de trouver une combinaison permettant 
de former un modéle couplé adéquat. Cependant, il ne s'agit pas d'un exercice trivial. 


Plusieurs approches sont présentées dans la littérature [51, 55, 56, 19]. 


2.4.2 Effets des modeles sur la précision des résultats 


Le choix du modele a un effet direct sur la précision des résultats. Par exemple en compa- 
rant des simulations en utilisant des modéles linéaires statiques, non linéaires statiques et 
non linéaires dynamiques, on peut obtenir des différences de l'ordre de 14-1696. Cet impact 
peut étre réduit si une stratégie non causale pour trouver l'optimum global est utilisée [6]. 
Le choix du modèle est donc une étape cruciale, puisqu'il affecte non seulement la précision 
des résultats, mais aussi le temps de calcul. Il faut donc trouver le bon compromis ou bien 


donner de la flexibilité à l'utilisateur. 


2.4.3 Validation du concept 


Une fois le concept choisi il faut le valider. De nos jours, il est pratiquement impensable 
d'aller directement à un prototype physique avant de l'avoir simulé. Il s'agit d'une étape 
cruciale avant la fabrication du bloc-batterie pour s'assurer qu'il respectera les limites 
d'opération et de sécurité nécessaire. Bien qu'il existe plusieurs modeles de bloc-batterie 
pouvant étre simulés comme décrit plus haut, une méthodologie avec une représentation 
plus graphique est de plus en plus populaire. La représentation énergétique macroscopique 
(REM) fut introduite il y a quelques années pour représenter graphiquement l'interaction 


de plusieurs sous-systémes reliés ensemble pour en simuler le comportement [57]. Cette 
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méthodologie permet principalement de modélisation les différentes interactions, mais aussi 
le modèle de contrôle ainsi que boucle de rétroaction. Cette méthodologie est un outil fort 


utile lorsqu'on doit modéliser des véhicules électriques [58]. 


2.5 Conclusion 


Le dimensionnement d’un bloc-batterie n’est pas une tâche simple. Beaucoup de variables 
doivent être considérées pour faire un choix éclairé. De plus, avec le nombre croissant 
de cellules disponibles, de comparer tout ce qui est disponible sur le marché avec toutes 
les combinaisons possibles devient une tâche impossible. À la suite de la revue de l'état 
de l’art, on peut remarquer plusieurs opportunités pour améliorer le dimensionnement de 
bloc-batterie. Tout d’abord, peu se sont concentrés sur la phase préliminaire de la concep- 
tion. L’obtention d’un ensemble de solutions faisables grâce à un problème de satisfaction 
de contraintes permettrait d’automatiser une partie du processus tout en analysant une 
multitude de combinaison possible. Ceci donnerait donc au concepteur une première vue 
d'ensemble sur ce qu'il serait possible de concevoir tout en pouvant obtenir une première 
comparaison de différentes cellules. Par la suite, on doit pouvoir distinguer une solution 
par rapport à une autre. Bien que plusieurs chercheurs se sont penchés sur l’optimisation 
du concept, il serait très intéressant de pouvoir évaluer la qualité de chacune des solutions 
faisables selon des critères définis. L'utilisation d'un problème d'optimisation contraint 
semble être la solution puisqu'il peut être considéré comme l'extension d'un probléme de 
satisfaction de contraintes, auquel on ajoute une couche d'optimisation. Dans l'optique 
oü l'on veut conserver le temps d'exécution rapide, le choix de la méthode d'optimisation 
est primordial. L'approche scalaire semble étre la meilleure solution pour obtenir les so- 
lutions non dominées dues à la nature discréte du probléme à résoudre, mais aussi pour 
conserver le temps d'exécution au minimum. Les 2 prochains chapitres présentent les deux 
parties majeures de ce projet de recherche. Tout d'abord dans le chapitre 3, le probléme de 
satisfaction de contraintes utilisé avec la représentation est présenté avec une utilisation 
dans un cas réel. Par la suite, le chapitre 4 présente la couche d'optimisation, le probléme 


d'optimisation contraint et la méthode de résolution avec démonstration sur un cas réel. 
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3.1.1 Resume 


Contenu : Cet article présente une nouvelle méthodologie de conception pour bloc- 
batterie de véhicule électrique et démontre les avantages avec un cas réel impliquant une 
moto électrique. Cette nouvelle approche propose de définir le bloc-batterie comme étant 
un problème satisfaction de contraintes (CSP). Au lieu d’itérer les concepts manuellement, 
le processus automatisé donne au concepteur la liberté d'explorer complètement l’espace 
de recherche. En intégrant le concept de contraintes dures et molles, cela aide à diriger 
la recherche dans la bonne direction. Surtout pour bien comprendre où le concept peut 
être flexible pour bien choisir le meilleur compromis. De plus, en priorisant les contraintes 
dans le bon ordre, le temps d’exécution peut être réduit de façon significative. Ce qui rend 
le processus encore plus rapide. Lorsqu’appliqués à un cas réel, les résultats démontrent 
que méme si le concept original était fonctionnel, des opportunités d'amélioration étaient 
toujours présentes. De plus, en comprenant quelles contraintes sont les plus restrictives, il 


aurait été plus facile de voir comment le concept aurait pu étre optimisé de facon efficace. 


Résultats : Les résultats démontrent que le probléme de satisfaction de contraintes (CSP) 
s'est avéré une méthode efficace pour obtenir l'ensemble des solutions faisables pour la 
conception d'un bloc-batterie. Également, lorsqu'appliqué à un cas réel, le concept obtenu 
démontre que des opportunités d'améliorations sont toujours présentes et permet de mieux 


comprendre ce qui contraint le concept. 


Contribution au mémoire : Les méthodologies de conception de bloc-batterie les plus 


répandus sont celles par essais/erreurs ou bien celles qui demandent des licences logiciels 
trés onéreuses. Il est donc trés difficile de trouver une méthode permettant surtout d'itérer 
les concepts en phase préliminaire pour comprendre l'effet des contraintes, mais aussi pour 
étudier les cellules disponibles sur le marché. De plus, l'implémentation d'une interface 
graphique facilite l'utilisation de cette méthodologie en plus d'avoir des temps d'exécution 
trés rapide. Finalement, l'utilisation d'un probléme de satisfaction de contraintes pour 
obtenir l'ensemble des solutions faisables servira de données d'entrée pour la méthodologie 


d'optimisation présentée dans le chapitre suivant. 
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This article presents a new framework for battery pack designs of electric vehicles and 
demonstrates the benefits with a real case study involving an electric motorcycle. The new 
approach proposes to define the battery pack as a constraint satisfaction problem (CSP). 
Instead of manually iterating the designs, the automated process gives the designer the 
freedom to explore the search space completely. By integrating the concept of hard and 
soft constraints, it helps to steer the search in the right direction. Mostly understanding 
where the design can be flexible to choose the best compromise. Moreover, by prioritizing 
the constraints in the right order, execution time can be reduced significantly. Making 
the process even faster. When applied to a real case scenario, results show that even if 
the original design was great, opportunities for improvement were still possible. Also, by 
understanding which constraints were the most important to prioritize, it would have been 


much easier to see where optimization in the design would have been the most effective. 


3.3 Introduction 


With the ever-growing electric vehicles (EVs) market, automotive constructors are now 
focusing on producing better, bigger and cheaper batteries. To be able to penetrate the 
mass market, they must produce batteries that can rival internal combustion engine per- 
formances. Range anxiety is probably the most common fear among the consumers along 
with charging time. Those factors must be considered and add complexity to an already 
non-trivial design. Many concurrent phenomena are occurring simultaneously making it 
very difficult to predict or emulate in computer simulations. There are also many param- 
eters to take in consideration such as cell chemistry, temperature of operation, capacity, 
power, etc. To optimize the design and to meet the requirements while maintaining low 
costs remains a complex challenge. Most of the research efforts are concentrated on the 
cell itself and not as much has been done regarding what happens when they are assem- 
bled together. The task of designing and assembling a battery pack is not trivial. With 
many interrelated variables involved, it becomes difficult to analyze this type of problem 
by hand. Nowadays, software and methodologies have been developed to facilitate this 
process. Companies like Altair®), ANSYS®), Gamma Technologies and others have devel- 
oped suites that allows the designer to build his own cells or battery packs but also allows 
to simulate them. They offer a lot of different features from design of experiments, opti- 
mization, simulation and data management [59], [4]. Their disadvantages are the prices, 
which are usually costly annual licenses plus the training and knowledge required to use 


the application. It is often required to have employees specially trained to use the product 
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to its full potential. In some cases, multiple software with multiple inputs and outputs 
are required to produce a design. Many authors have presented different methodologies to 
design battery packs. Over the last decade, optimization has become an ever-present solu- 
tion to improve one or more aspect of the design. Whether that is the thermal aspect [2], 
the cost |60], performance [61] or the mechanical aspect [5]. In [2], the authors present a 
multiobjective optimization for the thermal management of the battery. Using the particle 
swarm optimization algorithm, they were able to obtain the Pareto front considering the 
thermal behavior for the designs. The drawbacks are the simulation time which is around 
4 days and the need to possess a specific software license. A hybrid method is presented 
in [62], where the use of both gradient-free and gradient based method was utilized to 
accelerate the process and obtain better results. They optimized the design in respect 
to mass, volume and cost. The results obtained are promising, mostly for initial designs 
where subsequent analyses must be done to refine the solutions. But the complexity of 
the algorithm developed makes it very difficult to implement and use in industry. Another 
method is presented in [61], where models are used to optimize the design of the battery 
pack. Mostly focused on the thermal aspect, they are able to upgrade their initial design 
to a more robust design with lower temperatures. However, it is very difficult to analyze 
different configurations rapidly and is more focused on studying a particular configura- 
tion. Others have presented methodologies similar to these but integrating them within 
a tool. In [4], the authors present a tool that analyze the thermal behavior of the pack 
while studying different pack configurations. The tool allows more flexibility and ease to 
find the best configuration regarding the thermal aspect. The drawback is the complex 
architecture required to support the tool and the software license to run the computational 
fluid dynamics simulations. The mechanical aspect of the pack is often discarded in favor 
of thermal behavior. A methodology to optimize the layout and manufacturability of the 
battery pack is presented in |5]. It allows the designer to find the best mechanical layout 
and take in considerations some of the assembly process. But this method uses many hy- 
potheses and only takes in account the electrical aspect to model the pack. The thermal 
aspect is not considered in this approach. Many methodologies exist for the sizing of a 
battery pack and determining its configuration. They vary in complexity and exactitude. 
Plus, they all have a certain advantage over the others, but they are not perfect and op- 
portunities to better this process are still possible to explore. Some were highlighted in 
the last paragraph. This paper proposes a new framework to design battery packs based 
on constraint programming principles, by defining the problem as a constraint satisfaction 
problem (CSP) [63] and using finite domains constraints. This novel approach allows the 


designer to effectively and quickly explore the feasible set of solutions for the requirements 
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of the design. The paper is organized as follows: Section 3.4 introduces the constraint 
programming framework and explains how it can be used as a design tool for battery 
packs. Section 3.5 presents the case of study, the project specifications and a comparison 
of the results from the traditional method versus the new framework proposed. In section 
3.6, following the results shown in the previous section, the importance of constraints 
prioritization and the concept of soft and hard constraints is explained and their impact 


quantified. Finally, section 3.7 concludes with final remarks and future works perspective. 


3.4 Constraint Programming as a Design Tool 


Constraint programming is a programming paradigm where the interactions between the 
different variables are expressed in the form of constraints. The constraints that can be 
used may be of various kind, such as linear inequalities, non-linear or others. It has been 
used extensively in scheduling problems, but it is seeing more and more use in product 
design [64], [23]. This paper focuses on one aspect of constraint programming, named 
CSP where variables have finite domains and must respect several constraints. They are 
often of high complexity and require a heuristic or a defined search method to steer the 


exploration of the search space. Otherwise the computation time rises exponentially. 


3.4.1 Constraint Satisfaction Problem 


A constraint satisfaction problem is formally described by a triplet (X, D, C) which is 
defined as follows [23]: 


= X = (zi, £2, 23,..., Zn} is the finite set of variables of interest with n being the 
number of variables to be solved in the problem 

= D = (d, da, da, ..., v) is the finite set of domains for variables in X, where each x, 
of X gets its values in dn. 

— C = {c1, сә, C3, ..., Cm} is the finites set of constraints where m represents the number 


of constraints to be solved in the problem. 


3.4.2 Battery Pack Sizing Framework 


The battery pack design problem can be explained quite simply. There are variables 
which can take many values, but they must respect several constraints. The goal is then 
to find what are the possible combination of values for these variables where they will 
respect every constraint. Then, the search space can be minimized to only the possible 


solutions. This fits exactly the definition of a CSP. This concept is well known in the 
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artificial intelligence community. Very popular for solving scheduling problems, it has also 


seen some use in design problems and engineering [22]. 


The objective of the framework presented in this paper is to simplify the design process 
of battery packs. By modeling the pack design as a CSP, it facilitates and accelerates the 
process but also brings a flexibility not available with other methods. The designer enters 
his requirements, then the framework processes them in respect to the defined constraints 
and outputs the feasible set. Depending on the result, domains can be modified to be 
smaller or larger, variables can be changed, or constraints removed to explore other parts 


of the search space. The variables used in the framework are defined as follows: 


Design variables Symbol Unit Domain 
Number of cells in series N, cell [1, 250] 
Number of cells in parallel N, cell [1, 250] 
Distance between cells deell mm [1, 15] 
Cell model Cellmodey Brand and model LG HE2, LG HG2, 


LG MG1, LG MH1 
Table 3.1 Definitions of the Framework’s variables 


The specifications for each cell are defined in Table 3.2. 


Parameters Symbol LG HE2 LGHG2 LGMG1 LG MHI 
Diameter Dei 18.5 mm 18.5 mm 18.5 mm 18.5 mm 
Length Let 65 mm 65 mm 65 mm 65 mm 
Weight Meeu 44 g 44 g 44 g 44 g 
Minimum voltage Vmin,cell 2.5 V 2.5 V 2.5 V 2.5 V 
Minimum safe voltage Vis 3V 3V 3V 3V 
Nominal voltage лотва 9:00. V 3.66 V 3.66 V 3.66 V 
Maximum voltage Vmaz,cell) 4.2 V 4.2 V 4.2 V 4.2 V 
Capacity Qi 2.35 Ah 2.95 Ah 2.75 Ah 3.15 Ah 
Internal resistance Reell 80 mQ 80 mQ 120 mQ 120 mQ 
Maximum discharge current Imaz,cell 20A 20 A 10A 10A 


Table 3.2 Cells Specifications 


The constraints used in the framework are defined as follows: 

The constraint cı(Ns, Np, Ad.cı) is used to make sure that the design respects the max- 
imum volume available for the pack. Leeu, Deey and ‚Ad.eı are the length, the diameter 
and the distance between cells respectively. In this framework, all the available cells are 
Li-ion 18650 cells. Volmar is the maximum volume, in cubic meters, available for the 
battery pack and is defined by the designer. 
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Volmaz > N; d N, ` Гец H (Dea RE (2 E Adceu))? (3.1) 


The constraint ca(N,, Np) assures that the configuration respects the maximum mass al- 
lowed for the pack. Mu represents the mass of the cell, in kilograms, and depends on 
the selected cell model. Max is the maximum mass, in kilograms, allowed for the battery 


pack and is decided by the designer. 


Мааз > N, 2 Np t Меп (3.2) 
The constraint c3(Ns, Np) verifies that the configuration does not exceed the maximum 
price allowed for the design. Jeeu represents the price per cell and depends on the cell 


model and the market price. Jina, is the maximum price ($USD) for the battery pack 


and is determined by the designer. It only considers the price per cell. 


Jmaz > Ns t N, P ееп (3.3) 
The constraint c4(N,, Np) is used to verify that the configuration possesses enough energy 


for the usage of the pack. Үсеп is the nominal voltage of the cell and Emin is the 


minimum onboard energy, in kWh, necessary for the design and is decided by the designer. 


Emin < Ns d N, ` С сеп 3 Улот,сей (3.4) 


The constraint cs(Ns, Np) verifies that the configuration does not exceed the maximum 


allowed voltage Ving and where Vmaz,ceu is the maximum voltage of the cell. 


Vraz > М. Л Vmaz,cell (3.5) 


The constraint ce( N., Np) checks that the configuration has at least the minimum voltage 
necessary for the pack to operate. Where Vmin,ceu is the minimum voltage of the cell. For 


safety reasons, the cells are not allowed to go under 3V. 


Vinin < N; ў Vinin бей (3.6) 
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The constraint c7(Ns, Np) verifies that the pack can output the peak power required. 
Where Imaz,ceiı is the maximum discharge current of the cell and Рек the peak power 


needed for the application. 


Peak Ns E N, И Vnom;odit E Det (3.7) 


The framework also evaluates the maximum temperature that the battery pack will attain 
to give an idea of the thermal behavior to the designer. It is done with a simple steady 
state thermal estimation that considers the battery pack as a perfect cube, that the heat 
is uniform inside of the pack, only the convection is considered and is uniform on the 
six faces of the cube. The initial temperature is the ambient temperature. The model is 
defined by equations (3.8) to (3.10). Table 3.3 defines the values for the model’s variables. 


Variables Description Value Unit 
Rint batt Total resistance of the pack See 3.2 Q 
Tims Current representing the driving profile 58.3 A 
Tn Ambient temperature 25 °C 
L Length of the edge of the cube V Volume m 
A Surface area of the cube 6-1? m? 
Korick Thermal conductivity of PCM 7 ОЛ 
h Thermal convection coefficient 75 GR 


Table 3.3 Thermal Model Variables Values 


The heat generated during the driving cycle in the pack is calculated by (3.8). To be 
able to obtain the maximum temperature attained when in steady state, the total thermal 
resistance must be calculated with (3.9). Then, using the ambient temperature, the ther- 
mal resistance and the generated heat, the maximum temperature attained by the pack is 
obtained with (3.10). 


Q = Rint à I. (3.8) 

(2) 1 
Riot = 4 | 3.9 
© (Karick: A) (h- A) (99) 


docs — lamb t Riot А Q (3.10) 


3.5. CASE STUDY 39 


These constraints represent the common equations that need to be solved to start the 
design of a battery pack. By defining them as constraints, it becomes simple to change 
either the values of the requirements, the domains of the variables or the constraints 


themselves to change the search space. 


3.4.3 Framework Implementation 


The framework uses a generate and test approach to evaluate the CSP. Each variable of 
Table 3.1 will be iterated over its domain to evaluate each constraint. If all constraints 
are satisfied, the design is deemed feasible, otherwise it is discarded. These iterations 
will be done for each cell model and at the end, the thermal behavior will be calculated. 
When every configuration has been analyzed, graphs are drawn for each cell to be able 
to compare the differences. A graphical user interface has been created to facilitate the 


process. The implementation was done in MATLAB® for simplicity. 


3.5 Case Study 


To illustrate how this framework can be used in real applications, this section presents the 
analysis of a real case scenario. The battery pack studied in this paper was designed by 
the Electric Motorcyle Université de Sherbrooke (EMUS) team and can be seen in Figure 
3.1 It is a high-performance Li-ion battery pack built from the ground up for racing [19]. 


Figure 3.1 Université de Sherbrooke electric superbike prototype with high- 
lighted battery pack 


CHAPITRE 3. DESIGN D’UN BLOC BATTERIE HAUTE PERFORMENCE 
36 COMME UN PROBLEME DE SATISFACTION DE CONTRAINTES 


3.5.1 Project Specifications 


The goal was to design a battery pack for a racing application. For this case study, the 
powertrain and the chassis are already chosen. The superbike is powered by a 160kW 
EMRAX@®) 268 axial flux permanent magnet synchronous machine and is controlled by a 
Rinehart Motion Systems PM150. The chassis is an aluminium monocoque construction. 
These three components give the basis of what is required of the battery pack. By examin- 
ing the limiting values in the different datasheets and considering the racing environment, 


the requirements are formulated in Table 3.4. 


Requirements Value Unit 
Vols. 0.0694 m3 


Maar 120 kg 
dios 10000 $USD 
Emin 10 kWh 
Улат 700 V 
Vee 150 V 
Peak 160 kW 
Ps 80 kW 


Table 3.4 Battery Pack requirements 


3.5.2 Results and Discussion 


Now that the requirements are defined in Table 3.4, the framework can be used to analyze 
the different possibilities for the battery pack design. From Table 3.4, values are assigned 
to the different constraints. This enables the framework to evaluate the variables over 
their domains and to calculate the feasible sets. For the cell model, only the LG HE2, LG 
MG1, LG MH1 and LG HG2 were considered in this case study. Once every variable has 


been iterated over its domain, the Ns/Np graphs are produced. 


Figure 3.2 and Figure 3.3 are represented as number of cells in series on the Y axis and 
number of cell in parallel on the X axis. A configuration is represented as a Ns/Np couple. 
As stated in [12], the configuration choice is 135s 12p with LG HE2 18650 cells since 2016 
and was done manually. The first interesting thing is that Figure 3.2 shows that the design 
chosen can still be improved. For the same number of cells in series, 17 cells in parallel 
are required to respect every constraint of the project. Up from the original design of 12 
cells in parallel. But, if the LG HG2 is considered instead, the number of cells in parallel 
required is reduced because of the enhanced capacity of this model. Still keeping the same 
135 in series from the original design, 16 are now necessary in parallel. It is interesting 
to note that both the LG MG1 and LG MHI did not provide any feasible set. Mostly 


3.5. CASE STUDY 


37 


180 


160 


140 


2120 


100 


80 


60 


10 


15 


20 


25 


NP 


2 120- 


x 
HOO 


OO 


EEE 


IRRE 


15 


Np 


Figure 3.3 Feasible set of Ns/Np combinations for HG2 cell 
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because of their limitation on the maximum current and C-rate. Making them unsuitable 


for very high-power demand applications. In this case, the obvious choice is the LG HG2. 


It is of course a newer cell that was not available in 2016. If this fact is then considered, 


the LG HE2 was a good choice, but the configuration could still be improved. Instead of 
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the original 135s 12p configuration, the number of cells in parallel could be from 14 to 17 
depending on the number in series. To better fit in the operating voltage of the motor, 


the number of cells in series could be augmented significantly. 


When considering the mechanical packaging, adding more cells to the current design be- 
comes very difficult. Originally, the configuration of series and parallel cells was chosen 
for its ease of packaging. Adding more cells, either in series or parallel would be feasible 
but a difficult exercise. The framework provides a good indication of what is possible, but 
some work must be done afterward to obtain a real configuration. It allows the designer 


to see that some improvements are still possible. 
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Figure 3.4 LG HE2 thermal behavior of Ns/Np configurations 


Another important characteristic of the battery pack to consider is the thermal behavior. 
Ideally, the temperature that will be attained by the cells should be as low as possible. 
Both for safety and performance reasons. For the case studied in this paper, the method 
of cooling used is phase change material (PCM) from AllCell Technologies. The driving 
profile is extracted from data recorded at a real race that happened in 2017 at the New 


Jersey Motorsports Park circuit. Because of the simple static thermal model, the profile is 
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represented as the Irms seen by the battery. Figure 3.4 shows the maximum temperature 


of the battery pack for each of its possible configurations for the LG HE2 cell. 


By using the proposed approach, it is possible to obtain results that were otherwise very 
difficult to obtain manually. This gives freedom and flexibility in the design to facilitate 
the process and obtain better designs faster. This framework poses as special calculator 
helping the designer and does not replace him. Some insights and calculation are still 
needed to validate the design and transpose it to real life. But, as demonstrated in the 
case studied, it provides a better view of the space of possibilities and would probably 


have given a different direction in the manufacturing of the battery pack. 


3.6 Constraints Priorization and Relaxation 


In the previous section, the battery pack configuration analysis demonstrated the advan- 
tages of the new approach over a manual method. It is very difficult and time consuming 
to iterate over many designs possibilities by hand. Another thing that is difficult to con- 
sider when using a manual method is how to prioritize the different constraints and what 
are the limits of the search space. In a framework such as a CSP, constraints priorities can 
have a major impact on execution time. While evaluating a configuration, as soon as a 
constraint is not satisfied, the evaluation process stops and goes to the next configuration. 
By placing the most restrictive constraints first, execution time can be saved. To calcu- 
late every possibility for each variable over their domains takes around 835 milliseconds 
when using every constraint listed previously. When using only the two most restrictive 
constraints from Table 3.5, the calculation time goes up to 985 milliseconds. On the other 
hand, if the less restrictive constraints are used, such as cost, mass and minimum voltage, 
the total time of execution is 29,72 seconds. This shows the importance of choosing wisely 
how the problem is constrained and in what order they are evaluated. These calculations 
were made on an Intel®) Xeon®) E5-1650 3.50GHz with 32 Gb of RAM. 


The process of defining the priorities for each constraint will of course depend on the 
case studied. But most of the time, volume, power, energy, and voltage will be the four 
most restrictive constraints. For the case studied in this paper, the ranking is defined in 
Table 3.5. Once the constraints are ranked by priority, it is important to understand that 
they are not necessarily equally restrictive. Some will be hard constraint, where the value 
cannot in any case be exceeded. Some will be soft, where the value is more of a target 
that is malleable and has a lesser impact on the final design. Table 3.5 defines the type 


for each constraint for the case previously studied. 
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Rank Constraints Type 
1 Maximum volume Hard 
2 Minimum Peak Power Soft 
3 Minimum energy Hard 
4 Maximum voltage Hard 
5 Minimum voltage Hard 
6 Maximum cost Soft 


Table 3.5 Constraints Priorization 


It is important to understand how to effectively use these two types of constraints. Some- 
times, by relaxing the soft constraints, new designs become available that would otherwise 
have been missed. The importance of relaxation is primordial in these types of design prob- 
lems. With so many variables involved, the ability to relax some constraints to consider 


new possibilities becomes very useful in choosing the best compromise. 


Of course, the definitions of hard and soft constraints will be tied to the case studied. For 
some mass will be very important, for space vehicles for example, but in other case not so 
much. It then becomes important for the designer to understand these principles to define 
better requirements and constraints to be able to explore efficiently the search space of 


possible solutions. 


3.7 Conclusion 


In this paper, a new framework for battery pack design was presented. By defining the 
problem of designing a battery pack as a CSP, a systematic and formal approach was 
developed. This framework can be applied to any type of battery pack sizing. It allows 
the designer to define better requirements, but also better constraints for its design. It 
can help comprehend which constraints are the most limiting and then optimize around 
them to obtain better designs. The application of this framework on a real-world scenario 
shows that improvements opportunities can be detected as soon as at the beginning of the 
design process. It is also interesting to note that this framework could be applied to many 
other types of problem by adapting the constraints and parameters in short time. It is 


not limited only to high performance battery packs. 


For a more accurate exploration of the search space, new variables will be necessary. The 
form factor, of both the cell and the pack itself, will have a major impact on the final 
design and behavior of the pack. Sometimes using a prismatic cell instead of a cylindrical 
could save money or space and changing the form of the pack could have major impact on 
the thermal aspect. The depth of discharge could also be added to find the best window 
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for the application depending on the cycle applied to the pack. Finally, considering the 
packaging and the arrangement of the modules inside the pack is needed to better define 


which configurations are actually possible. 
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4.1.1 Resume 


Contenu : La conception de bloc-batterie pour véhicule électrique est un défi qui n’est 
pas simple à accomplir. Énormément de ressources informatiques, humaines et matérielles 
sont nécessaires. Beaucoup de temps et d’argent peuvent être économisés lors de la phase 
préliminaire du design en explorant le plus de solutions possibles dans le moins de temps 
possible. Cet article présente une nouvelle méthodologie de conception pour bloc-batterie 
qui permet des itérations rapides entre les concepts et de calculer les configurations op- 
timales pour donner le plus d'informations sur les meilleures solutions disponibles. Le 
bloc-batterie est représenté sous une configuration Ns/Np dans un probléme de satisfac- 
tion de contraintes (CSP) qui est résolu pour obtenir toutes les configurations possibles 
selon les requis spécifiques d'un véhicule. Le résultat est ensuite utilisé dans un probléme 
d'optimisation contraint (COP) pour déterminer les solutions non dominées selon les fonc- 
tions objectif sélectionnées. Les solutions non dominées sont obtenues en optimisant des 
sommes pondérées de fonctions scalaires pour des ensembles de poids dispersés, en plus de 
minimiser les métriques Li (distance Manhattan) et L (distance Chebyshev) par rapport 
à un point de référence, qui est établi avec des valeurs d'ingénierie. Ce point de référence 
permet de mettre en valeur les préférences du concepteur. Un cas réel est exploré oü cette 
méthodologie d'analyse est utilisée pour faire la conception d'un bloc-batterie pour une 
moto électrique. Il s'est avéré que cette méthodologie est trés efficace pour obtenir un 
prototype de bloc-batterie tout en comprenant quels sont les facteurs limitants ainsi que 
de quelle manière celui-ci pourrait être amélioré. Le bloc-batterie 245/36p conçu a ensuite 


été utilisé à la compétition MotoStudent 2018. 


Résultats : Les résultats démontrent que le choix d’une configuration séries/parallèles 
nest pas trivial. En effet, lorsque l’on considère le front de Pareto en entier avec des résul- 
tats bien dispersés, en projetant ce résultat sur une représentation Ns/Np, on comprend 
que de choisir le meilleur compromis n'est pas facile. Également, on démontre qu'une mé- 
thode haut-niveau, utilisant une combinaison de plusieurs sommes pondérées ainsi que des 
métriques Lı et Lo, permettent d'obtenir beaucoup d'informations sur le concept en peu 


de temps. 


Contribution au mémoire : Cet article présente la version définitive de la méthodo- 


logie de conception de bloc-batterie. En s'appuyant sur la représentation Ns/Np, l'ajout 


d'un probléme d'optimisation contraint ainsi que la méthode pour le résoudre, donne une 
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manière innovatrice d’obtenir de l’information sur le concept de bloc-batterie en phase 


preliminaire. 
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4.2 Abstract 


The design of electric vehicle batteries is a difficult challenge, requiring significant compu- 
tational, human and material resources. However, a lot of time and money can be saved 
in the early design phase by exploring an extended range of potential solutions in a small 
timeframe. This paper proposes a new framework to design a battery pack which allows 
quick iterations between designs and to compute optimized configurations to give more in- 
formation about the best available solutions. The battery pack is represented as a Ns/Np 
configuration in a Constraint Satisfaction Problem that is solved to exploit every possible 
configuration for specific vehicle requirements. The result is then used in a multi-objective 
Constrained Optimization Problem to determine the non-dominated solutions according 
to selected objectives functions. Non-dominated solutions are computed by optimizing 
weighted-sum scalar functions for well-dispersed sets of weights, as well as minimizing the 
Lı (Manhattan metric) and L4, (Chebyshev metric) distances to a reference point, which 
is established using practical engineering values. The use of weights and reference point 
allows to reflect the designer’s preferences. A case study is explored where this analysis 
framework is used to design a battery pack for an electric motorcycle. It proved to be a 
very efficient and effective methodology to obtain the first prototype of the battery whilst 
understanding the limiting factors but also how its design could be improved. The 24s/36p 
battery designed has been used successfully at the MotoStudent race in 2018. 


4.3 Introduction 


4.3.1 Motivation 


With the continuous growth of electric vehicles (EVs) sales in the automotive market, all 
the major companies are either developing an electric mobility strategy or starting to mass 
produce these vehicles. The main technology used in EVs is the lithium-ion battery, which 
has been evolving rapidly for the past two decades. As this technology is becoming more 
mature, the challenges are changing from how to connect elements together efficiently to 
make a battery pack to how can it be cheaper, more performant, smaller and lighter. 
While the technology for building a battery is now well known, this remains a non-trivial 
challenge. With several phenomena occurring simultaneously inside the battery, which are 
affected by the environment in which it operates, it is difficult to consider them all in the 
design phase. The optimization of every detail of the design takes a lot of computational, 
human and material resources to accomplish. Moreover, for every new design, each step 


must be repeated. This makes the process very difficult to be implemented in a fast, 
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efficient, cost effective and flexible way. Also, exploring the entire search space of feasible 
solutions becomes impossible for a human to do by trial and error. If the designer wants 
to evaluate every possible value for all the variables, the number of possibilities rapidly 
explodes. Therefore, this optimization process should be automated to reduce the time 
to market of battery packs. Several methodologies have been proposed in the literature 
and in the past few years companies have released commercial software suites to support 
that process. When compared to more traditional methods, these software packages clearly 
bring advantages but, to the knowledge of the authors, no methods proposed so far focuses 
on the early design phase to gain as much information as possible before going into costly 


optimization. 


4.3.2 Litterature review 


To properly size a battery pack is a tremendous challenge. It requires knowledge in many 
specific technical aspects such as mechanical, chemical and electrical engineering [19]. 
Moreover, many steps are required to go from requirements to a finished product. First 
deriving a number of good concepts during the early phase to then start optimizing the 
different aspects. Then, after many optimization passes, the first prototypes can be as- 
sembled and tested before reaching a product that will be used in a real vehicle. Several 
authors have presented different methodologies to size and design battery packs. Over the 
last decade, optimization has become an usual approach to improve one or more aspects 
of the design, including thermal [2, 65, 66], cost |60, 67], performance [61] or mechanical 
aspects [5]. In [2], the authors present a multi-objective optimization approach for the 
thermal management of the battery. The particle swarm optimization (PSO) algorithmwas 
used to obtain the Pareto front displaying the trade-off between minimizing the cell max- 
imum temperature, the power consumption of the thermal system and the area of one of 
the battery modules. The drawbacks were the simulation time, which is around 76 hours, 
and the need of a specific software license for COMSOL Multiphysics®), which can be 
costly. A hybrid method is presented in [62], where both gradient-free and gradient-based 
methods were utilized to accelerate the computation process and obtain better results. 
The design was optimized with respect to mass, volume and cost. The results obtained 
were promising, mostly for initial designs upon which subsequent analyses must be done 
to refine the solutions. However, the complexity of the algorithm developed makes it very 
difficult to be implemented and used in industry. Another method is presented in [61], 
where a electrochemical-thermal model was used to optimize the design of the battery 
pack. Mostly focused on the thermal aspect, the initial design was upgraded to a more 


robust design with a lower operating temperature in a race environment. However, this 
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methodology focuses on studying one particular configuration at a time, making it difficult 
to analyze different configurations rapidly. Other methodologies similar to these ones but 
implemented with customized software have been presented. In [4], the authors present a 
custom software implementation to analyze the thermal behavior for different pack config- 
urations. This software allowed more flexibility and eased finding the best configuration 
regarding the thermal evaluation aspect. The drawbacks were the complex architecture 
required to support the implementation and the need of a software license to run the com- 
putational fluid dynamics simulations. While volume, cost, mass or temperature are often 
considered as aspects to optimize, the mechanical aspect of the pack is often overlooked 
and thought about later in the process. A methodology to optimize the layout and man- 
ufacturability of the battery pack was presented in [5]. It allows the designer to find the 
best mechanical layout and take into consideration the assembly process. However, this 
method uses many hypotheses and only considers the electrical aspect to design the pack, 
the thermal aspect not being considered. Another interesting way of looking at sizing a 
battery pack is by sizing using the series-parallel configuration [21], i.e. displaying the re- 
sulting configuration in a simple way using the number of cells in series and parallel as the 
graph axes [19, 20]. This method of representation is becoming more popular, mostly in 
industrial settings for its easiness of understanding. As highlighted above, several method- 
ologies already exist for the sizing of battery packs. However, the development of a fast, 
flexible and reliable algorithm to tackle this problem poses many challenges. Multiple ways 
of representing the battery pack to be dealt with by an algorithm exist, which are mostly 
driven by the aspects that the designer wants to prioritize. The Constraint Satisfaction 
Problem (CSP) has regained a lot of popularity in recent years in engineering applications 
[22, 23]. Pioneered in the 1980s by artificial intelligence researchers [24], it is increasingly 
utilized to solve engineering problems. In CSP, given a set of variables each of them having 
a finite domain and interrelated by a set of constraints, the goal is to find feasible values of 
variables that satisfy every constraint. Several versions of a CSP exist. In the MAX-CSP, 
one seeks to satisfy the maximum number of constraints |25]. In the Weighted Constraint 
Satisfaction Problem (WCSP), one aims to minimize the total weight of the unsatisfied 
constraints |26]. While obtaining one or a particular set of solutions is of interest, optimiz- 
ing multiple conflicting objectives is very important to improve design options by enabling 
to consider the underlying trade-offs between distinct axes of evaluation. Because of the 
similarity of structure, the CSP definition can be extended to a Constrained Optimization 
Problem (COP) [68], which allows the consideration of one or more objective functions in 
the CSP results to find the best solutions. In general, multiple evaluations aspects that 


can be operationalized as objective functions are at stake to be optimized in a battery 
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pack. Different algorithmic approaches have been used to solve COPs enabling to cope 
with the combinatorial nature of the problem. Whenever mathematical programming al- 
gorithms cannot be applied due to the combinatorial dimension and/or non-linearity of 
optimization models, meta-heuristic approaches can be of help to provide good quality 
solutions with an acceptable computational effort; for instance, Tabu Search, working as 
a local search method, or Evolutionary Algorithms, enabling a global exploration of the 
search space by combining solutions with distinct characteristics [26]. While the goal of a 
COP is to find the optimal solution desirably with a lower computational effort, a good 
solution in a short time is often required for real-world applications [27]. In multiple ob- 
jective problems, a Pareto optimal front, or an approximation of it, should be computed. 
Multiple objective models and methods enable to consider multiple, conflicting and in- 
commensurate evaluation aspects of the merit of solutions seeking to obtain satisfactory 
compromise solutions belonging to the Pareto front, in the sense they establish acceptable 
trade-offs between the competing objectives. This front can be obtained computing one 
non-dominated (Pareto optimal) solution in each iteration, using scalarizing techniques 
such as weighted sums of the objective functions or minimizing a distance to a reference 
point the designer would like to attain, or by means of population-based approaches such 
as Evolutionary Algorithms in which a population of solutions expectedly converges to the 


true Pareto front (which, in general, is unknow in real-world problems). 


4.3.3 Contributions 


This paper proposes a new framework for lithium-ion battery pack sizing optimization 
using the Ns/Np representation presented ір [19, 20] in combination with a CSP approach, 
which is then extended to a multi-objective COP for which Pareto optimal solutions are 
computed by optimizing surrogate scalar functions. It is focused on the early stage of 
design to provide as much information as possible before moving to any expensive and 
specific optimization procedure. This proposed framework allows the designer to compare 
the results for different cells by providing series-parallel configuration graphs of those 
solutions. The thermal aspect is also considered, giving the designer further insights into 
the design. At the end, the framework outputs a set of non-dominated solutions evaluated 
with different metrics to obtain a broader picture of what suits the requirements best. 
Leaving space for the designer to choose according to its preference also. The design 
chosen can then be optimized further with more specific optimizations. The methodology 
presented in this paper is applicable to any type of batteries that are made of lithium- 
ion cells. As long as the driving cycle parameters and the battery packs requirements are 


known, any vehicle can be evaluated. To illustrate the potential practical usefulness of this 
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framework, a case study is presented to show how it can be used in a real case study where 


an electric superbike battery pack was built for an international student competition. 


4.3.4 Paper organization 


This paper is organized as follows. The motivation, literature survey, contributions, and 
organization have been addressed in this introduction. In section 2, the methodology is 
presented. A case study is presented in section 3. Results are presented and discussed in 


section 4. Conclusions are drawn and future work is outlined in section 5. 


4.4 Methodological approach 


The primary aims are to shorten time to market by reducing the time length of the first 
stage of the design process and obtain thorough results of the design before doing further 
engineering analysis that requires extensive resources and time. This section presents an 


overview of the methodology and details the different components. 


4.4.1 Constraint Satisfaction Problem 


There are multiple ways to define a CSP depending on the application. To explore all 
the possibilities for the design of a battery pack, some characteristics will inevitably be 
in conflict. This means that some trade-offs are necessary, which must be appraised and 
quantified. For these reasons, ће MAX-CSP definition is used [25]. This means that the 
algorithm will try to maximize the number of constraints satisfied and the quality of these 
solutions is measured by how many constraints, of equal importance, are satisfied. Using 
hard and soft constraints, the MAX-CSP is able to reject any non-feasible configurations 
and compute a score depending on the magnitude of the penalty associated with soft 
constraints that are violated. A hard constraint cannot be violated, and the corresponding 
design will not appear in the feasible set. A soft constraint may be violated, but a penalty 
score will be given depending on the violation degree and the design will still appear in 
the feasible set [69]. This framework implements a small modification of the MAX-CSP to 
feed the results directly into the COP algorithm to optimize the design. Since the designer 
already has some a a priori knowledge of the problem in the form of requirements from a 
client, this information enables to restrict the search space at the outset. The framework 
considers n variables, each of them having its own discrete finite domain. To solve the 
CSP, a generate and test approach is used, in which every configuration is calculated 


and evaluated. This ensures a complete evaluation of the search space to obtain as much 
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information as possible. For computation time reduction purpose, the number and domain 
range of variables should be the minimum required to fulfill the designer needs. This part 
of the framework is based on using the Ns/Np representation sizing methods presented in 


[8, 9], considering several improvements on constraint definitions, as detailed below. 


A constraint satisfaction problem is formally described by a triplet (X, D, C), which is 
defined as follows [24]: 
— X = {£1, 22, 23, ..., En} is the n-dimensional set of design variables of interest. 
= D = (di, ds, da, ...,2„} is the n-dimensional set of domains for variables in X, where 
each x; € X gets its values in a real - or set - valued d;, i = 1,...,n 
— C = {c1, со, C3,....,Cm} is the m-dimensional set of constraints defining the feasible 
region. 
For a solution to be considered valid to the CSP, it must be complete and consistent; 
i.e., a valid design solution must satisfy the defined domain of the variables and the hard 


constraints. 


The variables used in the proposed framework are defined in Table 4.1. 


Design variables Symbol Unit Domain 

Number of cells in series N, cell [1, 250] 

Number of cells in parallel N, cell [1, 250] 

Phase Change Composite ratio PCC atic % [10, 16] 
Cell model Cell,,44, Brand and model LG HE2, LG HG2, 


LG MG1, LG МНІ, 
Samsung 300 
Table 4.1 Definitions of the Framework’s variables 


The Phase Change Composite ! (PCC) ratio is set between 10% and 16% for design reasons 
and product availability. The cell model is considered a variable but is mainly used to 
define coefficients in the constraints. All cells listed in Table 1 are 18650 cylindrical cells. 


The specification parameters of each cell model are presented in Table 4.2. 
The constraints used in the framework are defined as follows. 


The constraints used in the framework are defined as follows: 
The constraint cı (Ns, Np, PCCratio) is used to make sure that the design respects the 
maximum volume available for the pack. The geometry and modularity of the pack is not 


considered in this constraint. This considers the volume of the cells plus the filling material 


1. PCC: Phase Change Composite thermal materials improve the safety and performance of lithium- 
ion battery packs. PCC is a lightweight combination of phase change material and graphite that absorbs 
and distributes heat [70]. 
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Parameters Symbol LG HE2 LGHG2 LGMG1 LG МНІ 30Q 
Diameter (mm) Den 18.5 18.5 18.5 18.5 18.5 

Length (mm) Lum 65 65 65 65 65 

Weight (g) Moen 44 44 44 44 45 
Minimum voltage (V) Е РЕТ 25 2.5 2:5 2:5 2:5 

Minimum safe voltage (V) Vente 3 3 3 3 3 
Nominal voltage (V) Vrom cei 3.66 3.66 3.66 3.66 3.66 
Maximum voltage (V) Vmanseei) 4.2 4.2 4.2 4.2 4.2 
Capacity (Ah) Ooh 2.35 2.95 2.75 315 29 

Internal resistance (mQ) Reel 80 80 120 120 85 
Maximum discharge current (A) Imaz,cell 20 20 10 10 15 


Table 4.2 Cells Specifications 


to handle as many cases as possible. The gap between the cell is filled with PCC from 
AllCell Technologies [70], which acts as the cooling medium. Leeu, Deen and PCC ratio are 
the length, the diameter and the ratio of PCC relative to the mass of the cells, respectively. 
In this framework, all the available cells are Li-ion 18650 cells. Volmar is the maximum 
volume, in liters, available for the battery pack, which is defined by the designer. To 
calculate the total volume of the pack, first the total volume of the cells is calculated by 
(4.1). Then, using the PCC ratio, the mass of the PCC, Masspco, is computed by (4.2) 
using the ratio relative to the cells mass. To obtain the volume of the PCC, its mass is 
divided by its density of 0.87 kg/L [71]. Finally, the total volume, in liters, of the pack is 
the sum of cells and PCC volume. The constraint c1(.N., Np, PC Cratio) is defined by (4.3). 


D 2 
Voleu = Ns- № - Гот ( =) (4.1) 
Masspcc = PCC Ratio x N, * N, * Meet (4.2) 
Masspcc 
Max = lee ( ) 4. 
Vol > Voleu + 0.87 ( 3) 


The constraintco((N,, Np, PCCratio) assures that the configuration respects the maximum 
mass allowed for the pack. Men represents the mass of the cell, in kilograms, and depends 
on the selected cell model. Mass,,,, is the maximum mass, in kilograms, allowed for the 
battery pack, which is defined by the designer. First the mass of the cells is calculated using 
(4.4). Using the result of (4.2), the total mass of the pack is computed by adding the cells 
total mass, Masseen, to the PCC mass Маѕѕрсс. The constraint co(Ns, Np, PCC ratio) is 
defined by (4.5), including (4.4) and (4.2). 
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Massi = N, E Np У Meel (4.4) 


Massmar > Masse + Mass pce (4.5) 


The constraint c3(Ns, Np) verifies that the configuration does not exceed the maximum 
price allowed for the design. Jeeu represents the price per cell and depends on the cell 
model and the market price. Jar is the maximum price ($USD) for the battery pack, 
which is determined by the designer. It is based only on the estimated price per cell. The 
constraint c3(.N,, Np) is represented by (4.6). 


doas = N; ү Np d dell (4.6) 


The constraint c4(N,, Np) ensures that the configuration possesses enough energy for the 
pack to be usable. Үсеп is the nominal voltage of the cell and Emin is the minimum 
on-board energy necessary, in kWh, which is established by the designer. Qeeu represents 
the capacity of the cell in mAh. The constraint c4(N,, Np) is defined by (4.7). 


Emin < Ns х № ; С) сеп Е Т (4.7) 
The constraint c;(N,, Np) verifies that the configuration does not exceed the maximum 


allowed voltage Vinaz, where Vmaz,ceu is the maximum voltage of the cell. The constraint 
cs(Ns, Np) is implemented by (4.8). 


Vraz > N; т Уа (4.8) 


The constraint ce( Ns, Np) requires that the minimum voltage of the pack at maximum 
power demand remains not lower than the minimum voltage required by the different 
components to function properly. Vince is the minimum voltage of the cell, Reen is 
the internal resistance of the cell and Гита is the maximum current demanded by the 
motor, which is divided by Np to obtain the current demand applied at the cell level. The 
constraint ce( Ns, Np) is defined by (4.9). 


Tee max 
Vmin < N; i (а = Reell t (Aem) (4.9) 
N, 
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The constraint c7(Ns, Np) verifies that the pack can output the peak power required. 
Imazx,cell is the maximum discharge current of the cell and Ppeak is the peak power needed 
for the application. The constraint cz(N,, Np) is stated by (4.10). 


Freak < Ns х № P ОРЕ, ї Tree (4.10) 


This framework uses a steady state thermal model to give an estimate of the maximum 
temperature reached by the pack during the cycle through resistive heating. The heat 


generation of a single cell qu is given by (4.11). 


cell == Reel $ 125; (4.11) 


To obtain Len , first the average power output of a cell is calculated with (4.12), which 
corresponds to the average power of the cycle Р, се average divided by the total number of 
cells. Then Fen is obtained with (4.13) by dividing the power output of a single cell Poen 


by the nominal voltage of the cell Үсеп: 


Роа 
pz _ 2 eyele,average 4.12 
ШИЛ? 5 
Pie 
Leu = ү £ (4.13) 
nom, cell 


Heat is generated nearly uniformly along the length of the rolled-up electrodes of the 
cell during discharge. The temperature gradient between the surface and the core is 
negligible, therefore allowing the lumped thermal model to be applied to the problem. 
This assumption applies to small diameters cylindrical cells, at moderate currents and 
under normal cooling conditions, which is the case in battery packs made of a multitude 
of cells enclosed in a cooling media. This allows the 3D multi-cell problem to be simplified 
to a one-dimensional single cell, surrounded by a layer of cooling material with steady heat 
losses. When assuming low contact resistance between the cell and the cooling media, both 
temperatures can be assumed identical. Thus, the thermal inertia of the pack mCp equals 
the sum of the thermal inertia of the cells Mep,ceıı and cooling media Mcp pee as given by 
equation (4.14). 


MCp = MCp, pce te Mep,cell (4.14) 
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À constant cooling power is assumed, may it be achieved via natural convection, forced 
air cooling or water cooling. It is possible to determine the temperature of the pack under 
steady-state conditions by computing the balance of heat generated by all of the cells беп 
and the amount that can be dissipated by the cooling system бе, over the cycle duration 
At. 


pack = cell Ne Np ка cool (4.15) 


Qpack == Qpack -At (4.16) 


Under those conditions, the temperature rise of the pack is obtained by dividing the total 
amount of heat to be stored over the thermal inertia, which is given by mo,. Using an 
ambient temperature Ty np) of 25 °C, the pack maximum temperature Tempmaz can be 
obtained with equation (4.18). 


= CO ack 


Mep 


Тетра = Tamo + AT (4.18) 


AT (4.17) 


4.4.2 Constrained Optmization Problem 


The main purpose of the application of the CSP is to remove all non-valid designs out 
of the results to allow for the consideration of optimized solutions which are feasible in 
the real world. To allow for the optimization of the previously obtained feasible set, the 
CSP framework is extended to a COP. Using the same mechanics as the CSP, the COP 
adds an optimization perspective, which is guided by the k-dimensional set of objective 
functions С = {91, 92, 93,- --, gk}, thus defining the quadruplet (X, D, C, G) enabling 
to handle different scenarios. The consideration of multiple objective functions allows for 
increased flexibility when optimizing the battery pack design by letting the designer define 
the objective functions in any combination to assess distinct scenarios. In this framework, 
six objective functions are considered, which can be chosen in any combination of three to 
facilitate the analysis. When specifying a given evaluation aspect as an objective function, 
it cannot be selected as a hard constraint. The definition and justification of the choice of 


the objective functions are presented in Table 4.3. 
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Obj; Objective function Definition Min or Max Justification 
Obj; Volmaz (4.3) Minimize Tight spaces 
Obja Мол (4.5) Minimize Lightweight vehicle 
Obj3 y орез (4.6) Minimize Cost 
Obj4 Erän (4.7) Maximize Autonomy 
Objs сай (4.10) Maximize Acceleration 
Obje Temps (4.18) Minimize Safety and aging 


Table 4.3 Definitions of the objective functions and the justification of their 
choice 


4.4.3 Solving the COP : Weighted-Sums, Lı and Læ Metrics 


Solutions to the multi-objective optimization problem are computed by means of solv- 
ing surrogate scalar problems with different technical parameters: weighted-sums with 
well-dispersed set of weights, as well as minimizing the Li (Manhattan metric) and Læ 
(Chebyshev metric) distances to reference points, the components of which are the levels 
in each objective function the designer would like to obtain. In addition to technical de- 
vices to introduce the designer's preference information into the computation process to 
identify suitable Pareto optimal solutions, the Loo metric enables to compute unsupported 
solutions (i.e., non-dominated solutions which are convexly dominated) besides supported 
non-dominated solutions resulting from optimizing weighted sum scalar functions and the 
Lı metric. Moreover, these different computation processes offer the designer different 
points of view of the possible trade-offs. Since some of the objective function values are 
dominant in terms of order of magnitude of their measurement scales, in comparison with 
other objective functions, a normalized value for each, normalizedyaiue,ooj, is obtained with 


(4.19) (for functions being minimized). 


тахітитоьу — valueos; 


(4.19) 


normalizedyalue,obj = - — 
MALİMUMopj — MINIMUM op; 


Equation (4.19) allows the objective function values to be normalized between 0 and 1. 
values; is the current value, minimum»; and maximum»; are the minimum and the 
maximum possible values of the considered objective in the feasible set. In order to offer 
the designer a complete overview of the non-dominated solution set, a regular grid of well- 
dispersed weights to be used in the optimization of the weighted-sum scalar functions has 


been defined (4.20). 
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k 
> = 1, where each w; varies in steps of 0.05 (4.20) 
j=l 


In (4.20), k is the number of objective functions selected by the designer and ш; is the 
weight value assigned to the corresponding objective function j. The idea is to obtain 
a more representative Pareto optimal front of the different trade-offs. The weights are 
positive to ensure Pareto optimality and they add-up to 1 for normalization purposes. To 
obtain the aggregate score of each design solution, the weighted-sum function (4.21) is used 
(it is considered that all objective functions are to be minimized, making the appropriate 
transformation in the maximizing objective functions, i.e., multiplying their coefficients 


by -1). 


min(F(x) = w- fi + we: fot... + we: fix) (4.21) 


s. t. : 4.3, 4.5, 4.6, 4.7, 4.10, 4.18 


The aim is solving problem (4.21) for each combination of weights to offer a comprehensive 


representation of the Pareto optimal front within the whole search space. 


Additional non-dominated solutions are computed to offer more information on the trade- 
offs in different regions of the Pareto optimal front, by computing the solutions that 
minimize the Lı and Læ distances to a reference point set by the designer. The scalar 
problem based on the L1 distance, the Manhattan metric, is defined in (4.22), where 
(21, 23; ..., 27) is the reference point for the design, which may be based on engineering 
values the designer would like to attain, k is the number of objectives to consider, and 
F;(x) is the objective function j. This is a compensatory metric since all the differences 
between each objective function value and the corresponding component of the reference 


point are considered. 


min( $ | zt — F(z) | ) (4.22) 


s. t. : 4.3, 4.5, 4.6, 4.7, 4.10, 4.18 


The scalar problem based on the Læ distance, the Chebyshev metric, aims to minimize the 


largest difference between the objective function values and the corresponding component 
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in the reference point is given by 4.23. This metric is non-compensatory because it only 


considers the largest deviation, i.e. the worst case. 


min( талуы, a(l — ©) |)) ЖҮ (4.28) 


s. t. : 4.3, 4.5, 4.6, 4.7, 4.10, 4.18 


This process allows an effective way of comparing different Ns/Np configurations, to un- 
derstand what the trade-offs between different designs are and how choosing different 
objectives can impact the final design. By using well-dispersed weights in weighted-sum 
scalar functions and minimizing the Lı and Loo distances to a reference point, a much 
broader and potentially useful Pareto front is generated, including the designer's require- 
ments and preferences (operationalized by means of weights and aspiration levels defining 
the reference point). In an industrial setting, this flexibility is a very important factor in 
making the early design process of battery packs more accurate and faster. The Pareto 
optimal front is also plotted on the Ns/Np graph to help in decision-making of which 
configurations is the most satisfactory for the designer. The main aim of this framework 
is not to give a final answer for the design of the pack, but rather offer the designer a flex- 
ible evaluation tool to assess different configurations and underlying trade-offs between 
the objective functions, which can then be arbitrated using the designer's preferences or 


further requirements. An overview of the framework algorithm is presented in Figure 4.1. 


4.4. METHODOLOGICAL APPROACH 


99 


Requirements 


CSP 4 


Display feasible set for each selected 
cell 


point 


Hard constraints and optimization 
objectives selection 


objective 
elected? 


Driving cycle average 
power and duration 


Normalize each objective value for 
every configurations in the feasible set 


Calculate the Manhattan and 
Chebyshev distances to the reference 
point for each configurations in the 
feasible set 


Weighted sum calculation for each 
configurations with scattered weights 


Find minimum value of the weighted 
sum for each set of weights and for 
the Manhattan and Chebyshev 
distances 


Output the non-dominated solutions 
onto the Ns/Np plot with the optimal 
solution highlighted. 


Constraints 
relaxation 


ctory? 


Selection of final configuration 


\ 
\ 
\ 
\ 
\ 
N 
М 


N YES 


\ 
\ 


А EN 
Is configuration. 

N satisfa 
\ 


PO 


Figure 4.1 Design framework flowchart 


COP 


A Graphical Designer Interface (GUI) is developed in MATLAB to facilitate the interac- 


tion. Its primary goal is to configure the framework to the designer's preferences regarding 


requirements and optimization objectives. The GUI also allows the designer to save and 


load configurations, use a database of cells for comparison, visualize Ns/Np graph of the 


feasible set and the temperature and optimize his designs. The final interface can be seen 


in Figure 4.2. 


The next section presents the case study where the framework has been used to design a 


battery pack for an electric motorcycle. 
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Figure 4.2 Framework GUI to design multi-cell battery pack 


4.5 Case Study - Electric Motorcycle Prototype 


To illustrate how this framework works for real applications, the battery pack of the 2018 
EMUS MotoStudent electric race bike is used in this section. Each step used by the 
framework will be demonstrated to understand how the framework can help the designer 
in the decision-making process. For the MotoStudent race, the competition organizers 
give every team a permanent magnet synchronous motor operating at a maximum voltage 
of 110 V, which outputs 48 kW of peak power at 600 A and it is controlled by a SEVCON 
Gen4 Size6 motor drive. For ease of assembly, modularity, reliability and safety reasons, 
the pack must be divided into 6 modules. The total number of cells must be a multiple of 
6. The custom chassis is an aluminum monocoque construction, allowing maximum space 
for the battery pack. These components are fixed and therefore, for this case study, the 
motor, the controller and the chassis cannot be modified. The goal is then to define the 
requirements for the framework to create the battery pack and optimize its configuration 
based on these components, as shown in Table 4.4. T'hese requirements are used to define 


the constraints and their limit values. 


To help in evaluating and comparing the different possible Ns/Np configurations, a driv- 
ing cycle was simulated to represent a possible MotoStudent race. For this cycle, when 
excluding the high transients demand of power, the estimated average continuous power 


is around 23 kW over 700 s. These two parameters are used by the thermal model to 
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Figure 4.3 The 2018 EMUS electric race bike on the Race Track at Circuit 
Motorland Aragon, Spain. Credits for Photo : MotoStudent Organization, Jose 
Espallargas and asm press, 2018 


estimate the maximum temperature attained by the pack during the cycle. A part of the 


driving cycle is illustrated in Figure 4.4 and is obtained by the method used in [19]. 


The framework is configured for the case study of this battery pack as presented in Table 
4.4. 


Constraints Value Unit Hard constraints Objectives 


Volmaz 40 L Yes No 
Mmazx 60 kg No Yes 
dius 5000 $USD No Yes 
Emin 9 kWh Yes No 
Vic 110 V Ves No 
Vmin 34 V Yes No 
Ppeak 42 kW Yes No 
Pont 20 kW Yes No 

T'empgas 60 °C No Yes 


Table 4.4 Configuration for the framework selected to solve the problem 
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Figure 4.4 Driving cycle for the MotoStudent race representing ап estimated 
typical lap at Circuit Motorland Aragon 


The reference point used to calculate the solutions that minimize Manhattan and Cheby- 
shev distances is defined as follows. For each objective, a reference value is chosen accord- 
ing to an ideal value for the design of the battery pack, which may be very ambitious in 
real life, based on engineering calculations. For the mass, the value is based on the energy 
density of the Samsung 30Q cell and the energy required for the race. In this case, the 
cell has an energy density of 235 [71] and with the energy requirement of 9 kWh gives 
a mass of about 38 kg as the reference value. The maximum temperature is based on the 
optimal temperature of operation of the cell as per the cell datasheet, which is 25°C [20]. 
Finally, the cost is the budget for the battery pack that the team had at the time of the 
design to buy all cells, which is set at 4300$ USD. 


4.5.1 Results and Discussion 


In the case study described in the previous section, the battery pack design is assessed 
using the evaluation aspects selected in Table 4.4, to minimize mass, temperature and 
cost while considering a racing application. To demonstrate how the proposed framework 
can be used, first the feasible set graph is presented which represents only the results of 
the CSP dysplaying how constraints play a role in defining the feasible set. Then a step 
by step analysis of the optimization process starting with one objective and going up to 
three objectives is performed. While this process is normally executed at once, the results 


can be understood more clearly by presenting how the framework works with the data. 
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Although results for other cells have been obtained, only the output for the Samsung 30Q 


is shown in this section for the sake of keeping the paper concise. 


4.5.2 Feasible set of Ns/Np Configurations 


Once the framework is parameterized, the first step is to solve the CSP to obtain the 
feasible set of Ns/Np configurations for each one of the chosen cells. Figure 5 shows the 
feasible set obtained for the Samsung 30Q. This represents each series/parallel configu- 
ration that does not violate the hard constraints selected in Table 5 for the design and 
are plotted using the Ns/Np representation. To understand the shape of the feasible set, 
Figure 5 also displays how each hard constraint is limiting the search space. The solid 
lines represent different limit values that the solution must respect. Although there are 
more than 4 constraints used in this case, only the most restricting ones are shown in 
Figure 4.5.Once all the non-feasible designs are removed and many possibilities are still 
viable, the next step is to solve the COP for one objective using the CSP results as input, 


in this case minimizing mass. 


4.5.3 One Objective - Minimizing Mass 


Once the feasible set is computed, it can be fed to the optimization process. It serves as 
the starting point for the COP, providing the configurations that respect the requirements 
for the vehicle. For illustration purposes, each step will be evaluated separately. When 
optimizing for only one objective, there is a single optimal solution (although alternative 
optima may exists in some cases). The configuration with the lowest mass is 24s/35p. 
Since mass is mostly dependent on the number of cells, the configurations with the lowest 
masses will be those that first minimizes the number of cells in parallel and then in series. 
In the next subsection, two objectives are evaluated in the COP and the changes in the 


results are highlighted. 


4.5.4 Two Objectives - Minimizing Mass and Temperature 


In Figure 4.7 a), the Pareto front graph shows a more interesting picture when comparing 
two conflicting objectives. By adding a new objective, the algorithm tries to find a bal- 
anced solution between minimizing the mass, and the maximum temperature of operation. 
By using well-dispersed weights for the optimization of weighted-sum scalar functions, a 
broader set of design options is offered as displayed by the Pareto optimal front. The 
solutions obtained by minimizing the Lı and L,, metrics bring more information thus 


offering further insights for the decision-making process. For this purpose, the reference 
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Figure 4.5 Series-Parallel configurations feasible set and a representation of all 
the constraints and how their limits affect solving the CSP. 


values for the maximum temperature and mass should be established. Looking at Figure 
4.7 b), the choice is not trivial. Considering the Manhattan distance, the best solution 
is 25s/40p, while for the Chebyshev distance is 24s/38p; other non-dominated solutions 
highlighted in red are still good candidates embodying different trade-offs between the 
two objective functions. The 24s/38p design configuration would be a good choice for 
modularity reasons. For the race bike, the pack is divided into 6 modules to form the 
whole battery, thus meaning that the total number of cells must be a multiple of 6 to fit. 


Each color represents the same configuration in both figures. 


To better understand the results of the temperature optimization, Figure 4.8 presents the 
maximum temperature reached by each possible configuration. The configuration 24s/38p 
has a maximum temperature of 51 °C. This is well below the safety limit that is set at 


60 °C, which means there could be room to reduce mass if it is acceptable to have a 
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Figure 4.6 Feasible set with results of optimization, configuration with lowest 
mass highlighted 24s/35p in red 


higher temperature [72]. This allows to match the results of Figure 4.7 b) with the Ns/Np 
representation in Figure 4.8. This can be used posteriori to understand how to obtain 
better results in a subsequent iteration, if the temperature is not suitable. Considering 
an industrial setting, these graphs allow for better communication and understanding 
of design choices. The next subsection presents the evaluation with 3 objectives: mass, 


temperature and cost. 


4.5.5 Three Objectives - Minimizing Mass, Temperature and Cost 


The same process is repeated for three objectives, which are minimizing the mass, minimiz- 


ing the maximum temperature of operation and minimizing the cost. The cost objective 
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goes in the same direction of the mass objective, i.e. to reduce the number of cells. In 
Figure 4.9 a), the shape of 3D Pareto set is very similar to the one obtained with two 
objectives as shown in 4.7 a). This is mostly due to the relation between cost and mass, 
being both dependent on the number of total cells. It also illustrates that to be able to 
lower the maximum temperature of the pack, both mass and cost will increase. The results 
of the optimization process offer novel insights. Modularity becomes less of an issue for 
the solutions obtained by minimizing the L4, and Læ distances because they consist of 24 
cells in series. The decision for the final configuration to be chosen should also be based on 
what a racing electric race bike needs to perform. For performance purposes, the electric 
bike must be as light as possible; thus, the 24s/36p design solution could be the selected 


configuration for this application. 
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Figure 4.7 Results obtained from solving the COP with 2 objectives (Mass 
vs. Temperature): (a) Pareto optimal front (b) Feasible set with optimization 
results highlighted. 


In Table 4.5, a summary of the 24s/36p design configuration is presented. 


It is interesting to note that all this information can be obtained very quickly with this 
framework, thus enabling to assess different scenarios in terms of input coefficients and 
preference parameters (e.g., the components of the reference point). When evaluating 6 
different cells, the resulting computing time is 6.55s, which means that this framework 
can be integrated in an industrial workflow easily. This calculation was made on an 
Intel® Core™ i7-6500 CPU at 2.50GHz with 8 Go of RAM laptop. Even if iterations 
can be done quickly, it is important to approach the design with intent and gather the 
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Figure 4.8 Maximum temperature of each Ns/Np configurations based on the 
driving cycle parameters 


relevant information. Each objective selected has its own behavior and shifts the set 
of non-dominated solutions according to its own optimization direction. Therefore, the 
trade-offs between the objectives considered should be carefully assessed to obtain the 
best design results. The proposed framework is designed to handle a flexible number of 
objectives combination and it is to the designer’s discretion to select the more adequate 


combination. 
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4.6 Conclusion 


This paper proposes a new framework to optimize the sizing of battery pack design in the 
early phases of the design process. Using a combination of CSP/COP, multiple scalarizing 
weighted-sums, Lı апа La distance functions to compute non-dominated solutions to a 
multi-objective optimization problem, this approach enables the presentation of different 
configurations in a very visual and designer-friendly manner. The CSP allows a fast and 
efficient way of removing Ns/Np configurations that would not be possible or interesting 
to build. Then, the COP allows us to evaluate the quality of each possible configurations 
by highlighting the Pareto optimal front using different computation techniques to provide 
a comprehensive view on the possible trade-offs. Moreover, the Ns/Np temperature graph 
is a novel way to assess the different configurations. In just a glance, the difference 
between the configurations in terms of temperature can be evaluated without the need 
of experiments or complex models. Using the configuration specified in section 3, the 
computing time is 6.55s to evaluate 6 cels. As shown in the case study, the framework 
was used to design a battery pack for an electric motorcycle. This framework requires no 


costly software licenses and can be used on any computer. In addition, the way the Pareto 
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Figure 4.9 Results obtained from solving the COP with 3 objectives 
(Mass/Temperature/Cost): (a) Pareto optimal front for the 3 objectives (b) 
Feasible set with optimization results highlighting the solution 24s/36p 


optimal front is constructed and highlighted in the Ns/Np graph shows that choosing 
a configuration is not a trivial endeavor. Future works will add more variables to the 
framework and a revamped GUI. Also, packaging should be considered when calculating 
weight, volume and cost to give a more realistic design. Being able to consider how the 
pack is divided in modules in such a framework is also an interesting subject that needs 
to be researched. Finally, comparing the results obtained with the proposed framework 
with the ones from a genetic algorithm for problems with a higher number of variables 
(which may lead to combinatorial difficulties) is an interesting methodological topic for 


future works. 
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Parameter Value Unit 
Number of cells in series 24 - 
Number of cells in parallel 36 - 
Total number of cells 864 - 
Volume 34.87 L 
Mass 62.8 kg 
Cost 4320 $USD 
Energy 9.64 kWh 
Maximum voltage 100 V 
Minimum voltage 42,8 V 
Peak power 46 kW 
Maximum temperature 57.85 °C 


Table 4.5 Summary of final Ns/Np configuration chosen for the battery pack 
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CHAPITRE 5 


Conclusion 


5.1 Sommaire des travaux 


Les travaux présentés dans ce mémoire portent principalement sur deux aspects du dimen- 
sionnement de bloc-batterie lithium-ion : l'automatisation d'une partie du processus ainsi 
que l'évaluation de la qualité des solutions non dominées. La combinaison d'un probléme 
de satisfaction de contraintes et d'un probléme d'optimisation contraint a permis de réunir 
l'obtention des solutions faisables avec une méthode d'optimisation pour obtenir des ré- 
sultats concluant extrémement rapidement. De plus, la méthode d'optimisation employée 
a permis de conserver les temps d'exécution sous les 10 secondes tout en obtenant un front 


de Pareto contenant des solutions non dominées plus diversifiées. 


Le chapitre 3 présente la premiére partie de la méthodologie soit le probléme de satisfaction 
de contraintes. Ceci a permis de démontrer qu'il existe une maniére simple et efficace de 
retirer toutes les solutions ne correspondant pas au requis et gráce à la représentation 
série-paralléle, le concepteur obtient un graphique permettant d'analyser rapidement les 
configurations possibles. Cette manière innovatrice de présenter le probléme de conception 
de bloc-batterie permet de facilement obtenir l'ensemble des solutions faisables avec un 
faible coüt en ressource informatique. Cette méthode peut étre facilement modifiée pour 
inclure de nouvelles variables, contraintes ou objectifs pour bonifier celle-ci. Elle permet 
aussi de comparer différentes cellules et obtenir des indications sur le meilleur candidat 
potentiel selon les requis. Pour démontrer le potentiel de cette méthode, elle a été utilisée 


sur un cas réel de dimensionnement de bloc-batterie pour une moto électrique. 


Par la suite, les travaux décrits au chapitre 4 ont démontré qu’il était possible d'opti- 
miser le concept de bloc-batterie en phase préliminaire suffisamment pour avoir une idée 
claire du concept potentiel pour les prochaines phases. En effet, l'utilisation du probléme 
d'optimisation contraint se servant du résultat du probléme de satisfaction de contraintes 
comme données d'entrée a permis d'utiliser les résultats du chapitre 3 comme base pour 
finaliser la méthodologie. De plus, en utilisant plusieurs sommes pondérées de fonctions 
scalaires ainsi que les métriques Lı et Læ, ceci a permis l'obtention d'un front de Pa- 
reto large et diversifié. Pour pouvoir gérer les différentes unités et les fortes tendances de 


certaines variables, la normalisation des valeurs des objectifs s'est prouvée trés efficace. 
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Il s’agit d’une manière novatrice d’optimiser le dimensionnement d’un bloc-batterie. Bien 
que la majorité des méthodes présentées dans l’état de l’art se concentre sur un concept 
en particulier ou tente d'optimiser entre un et trois paramètres, aucun ne se concentre 
sur la phase préliminaire. En passant d'un probléme multi-objectifs à un probléme mono- 
objectif, ceci permet de garder le temps d'exécution rapide sans augmenter la complexité. 
La conjonction de plusieurs sommes pondérées avec des poids variés et les métriques on 
permit de quantifier la qualité des différentes solutions non dominées obtenues. Puisque 
ces métriques utilisent un point de référence basé sur les préférences de l'utilisateur, ceci 
permet d'intégrer les préférences de celui-ci facilement tout en orientant les résultats dans 
une direction plus favorable pour lui. Un cas de dimensionnement réel a été effectué avec 


cette méthodologie pour illustrer l'efficacité et le potentiel de cette nouvelle méthodologie. 


5.2 Travaux futurs 


Le projet de recherche décrit dans ce mémoire permet l'ouverture de plusieurs voies de re- 
cherches et d'améliorations possibles. Tout d'abord, déterminer une maniére de considérer 
la structure du bloc-batterie dans le calcul de volume, de masse et de coüt pour obtenir 
un résultat plus prés de la réalité. Bien que déterminé la configuration de la batterie soit 
une étape cruciale, il ne faut pas oublier que la structure entourant celle-ci est un défi 
considérable à entreprendre. De pouvoir considérer ce facteur dés la phase préliminaire 
permettrait d'obtenir un portrait complet du produit final sans avoir à le construire ou le 
modéliser dans un logiciel. Il serait intéressant de faire le méme exercice pour inclure la 


masse, le volume et le coüt des composantes électroniques. 


Le dimensionnement ne se limite pas au véhicule. L'utilisation de cette méthodologie sur 
des applications stationnaires pourrait amener des résultats trés intéressants. En compa- 
rant avec des résultats obtenus avec les méthodes traditionnelles, on pourrait confirmer 
que cette méthodologie peut s'étendre à tout type de bloc-batterie, si les requis sont définis 


selon les intrants demandés. 


Bien que l'interface graphique facilite grandement l'utilisation de cette méthodologie, le 
développement d'une application compléte pourrait améliorer grandement la convivialité 
de celui-ci. Par exemple, en ajoutant une base de données contenant les données de dif- 
férentes cellules de différents manufacturiers. De plus, d'autres variables, contraintes et 


objectifs pourraient étre ajoutés. 


Finalement, il serait intéressant de comparer les résultats obtenus avec la méthodologie 


présentés avec un algorithme génétique en intégrant un plus grand nombre de variables. 
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La méthode exhaustive possède des limites en termes d’exécution lorsque le nombre de 
variables augmente de façon significative, l’utilisation d’un algorithme pour passer outre 


cet obstacle serait clairement utile dans une version future. 


76 


CHAPITRE 5. CONCLUSION 


ANNEXE A 


Guide d'utilisation de l'interface graphique 


A.1 Procédure d'ajout d'une nouvelle cellule 


L'interface graphique permet au concepteur beaucoup de flexibilité. Que ce soit au niveau 
des paramétres évaluer que les cellules disponibles. Lorsque celui-ci regoit une nouvelle 
cellule et la teste pour ensuite vouloir l'utiliser dans de futurs designs, il peut l'ajouter 
directement dans l'interface en suivant la procédure suivante. Un exemple utilisant la cel- 
lule NCR18650PF de Panasonic est démontré ci-dessous dans l'environnement MATLAB. 
Tout d'abord, les caractéristiques suivantes sont extraites de la fiche technique. 


Variable Valeur Unité 
Nom NCR18650PF string 
Capacité 2,1 Ah 
Courant maximal 10 A 
Crate 3,7 C 
Résistance interne 0,043 Ohms 
Tension minimale 2,5 V 
Tension nominale 3,6 V 
Tension maximale 4,2 V 
Diamètre 0,0185 décimétres 
Longueur 0,065 décimétres 
Masse 46,5 grammes 
Coüt approximatif 4 $USD 


Tableau A.1 Caractéristique d'une cellule lithium-ion 


Idéalement, certaines données comme la capacité, la résistance interne et la masse peuvent 
être obtenues en laboratoire pour une précision accrue. À noter que les unités sont impor- 
tantes puisque l’algorithme se base sur celles-ci pour faire les calculs. Ces données sont 
ajoutées au fichier loadCell.m comme présenté à la Figure A.1. 


La prochaine étape est d'ouvrir le fichier UI_v_0_1.fig avec GUIDE dans MATLAB. П 


faut éditer la liste de cellules dans le bas à gauche de la Figure A2. 
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Editor - CA 


\louis\ PDT_1 k_Sizing\loadCell.m* 


EDITOR PUBLISH VIEW 
+ J a LG] Find Files Ss men Б) f E - [> ie ТАРИМ > 
New Open Save ME ARMEE ү Breakpoints Run  Runand |=} Advance Runand 
- ох (XP ov Find ~ тает [5] $3 + М ~ Advance Time 
FILE NAVIGATE EDIT BREAKPOINTS RUN — — | 
DistanceOptimizationProcess3Obj.m | DistanceOptimizationProcess2Obj.m | ThermalOptimizationBrute.m Ul v. O 1.m WeightedSumOpt3Obj.m loadCell.m* | CSI 
E RECONNU VU UT LRU " И " 
Se CurrentCell.R theta-1./0.018; 
[iro CurrentCell.Max Current - 26; 
Lire case 8 %SAMSUNG 30Q 
HE CurrentCell.Name = "Samsung 300'; 
119 — CurrentCell.C cell = 3.09;%[Ah] 
120 — Currentcell.Max_Current = 15; 
BEZ CurrentCell.Crate- CurrentCell.C cell / CurrentCell.Max Current; % [С] 
ph CurrentCell.R int- 0.085;%[Ohm] 
123 — CurrentCell.Vcell min-2.5;$[V] 
124 — CurrentCell.Vcell nom-3.61;$[V] 
255 CurrentCell.Vcell_max=4.2;%[V] 
i26 — CurrentCell.D cell-18.5e-2;$[dm] 
RIS CurrentCell.l cell-65.2e-2;$[dm] 
128 — CurrentCell.m cell-0.045;$[kg] 
129- Currentcell.J_cell=5;%[$] 
ізо — CurrentCell.wt-2e-3;$ [m] 
Bez CurrentCell.m pcm-0.0075;$[kg] 
Bae CurrentCell.R_theta=1./0.018; 
B3- case 9 %NCR18650PF 
134 — CurrentCell.Name = 'NCR18650PF'; 
fasi CurrentCell.C cell = 2.7;%[Ah] 
136 = CurrentCell.Max Current - 10; 
BIS CurrentCell.Crate- CurrentCell.C cell / CurrentCell.Max Current;$[C] 
138 — CurrentCell.R int- 0.043; %[Ohm] 
lg CurrentCell.Vcell min-2.5;$[V] 
L40 — Currentcell.Vcell_nom=3.6;%[V] 
141 — CurrentCell.Vcell max-4.2;$[V] 
Malz CurrentCell.D cell-18.5e-2;$[dm] 
АЗЕ CurrentCell.l се11=65.2е-2; % [dm] 
144 — CurrentCell.m cell-0.0465;$[kg] 
145 — CurrentCell.J се11=4;+[5] 
146 — otherwise 
147 — CurrentCell.C cell-2.85;$[Ah] 
148 — СиггепЕСе11.Сгаїе=10; % [С] 
149 — CurrentCell.R int-0.02;$[Ohm] 
150 = Currentcell.Vcell_min=2.5;%[V] 
Bez CurrentCell.Vcell nom-3.7;$[V] 
525 CurrentCell.Vcell max-4.2;$[V] 
53 = CurrentCell.D се11=18.5е-2; % [ат] 
[Ба = CurrentCell.l cell-65.2e-2;$[dm] 
655 Currentcell.m_cell=0.044;%[kg] 
156 — CurrentCell.J cell-2.25;$[$] 


Figure A.1 Fichier loadCell.m 


Ensuite, dans l'éditeur de propriétés, ajouter le nom de la cellule au bas de la liste comme 
à la Figure A.3. 


Sauvegarder les modifications et relancer l'interface graphique. La nouvelle cellule sera 
maintenant disponible pour étre utilisée. 


A.2. SAUVEGARDER ET CHARGER UNE CONFIGURATION 
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Figure A.3 Éditeur de propriétés de la liste de nom des cellules 


A.2 Sauvegarder et charger une configuration 


Si le concepteur a effectué un dimensionnement à l’aide de l'interface et veut conserver la 
configuration de celle-ci, il peut le faire avec l’option Save Scenario. Celle-ci se retrouve 
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dans le menu Save en haut à gauche de l'interface. Une boite de dialogue ouvrira automa- 
tiquement et la configuration pourra étre sauvegardée à l'endroit voulu sur l'ordinateur. 


Lorsque le concepteur voudra charger cette configuration, il peut le faire avec l'option Load 
Scenario dans le menu Load en haut à gauche de l'interface. Une boite de dialogue ouvrira 
automatiquement pour sélectionner le fichier de configuration voulu. Celui-ci configurera 
automatiquement l'interface. 
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